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^ . ABSTRACT 



An introduction to modern cosmology is given, focusing on theories for the 
origin of structure in the Universe. After a brief review of the theory of growth 
of cosmological perturbations and a summary of some important observational 
results, most of the article is devoted to discussions of the inflationary Universe 
scenario and of topological defect models of structure formation. The final chapter 
is a summary of the theory of cosmic microwave background temperature fluctua- 
tions and of the current observational status. 
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1. Introduction 

Cosmology has over the past fifteen years emerged as a vibrant and exciting 
subfield of physics. It is based on the marriage of quantum field theory and particle 
physics on the one hand with classical general relativity on the other. One of the 
main goals of modern cosmology is to explain the structure of the Universe on 
the scale of galaxies and beyond. Thus, the experimental/observational basis of 
the field lies in astronomy, and there is a lot of interaction between theoretical 
cosmologists and observational astronomers and astrophysicists. 

The main goal of these lectures is to give an introduction to the two most 
developed classes of structure formation theories: those based on infiation and those 
based on topological defects. I will give a brief survey of relevant observational 
results from large-scale structure surveys and from searches for cosmic microwave 
background (CMB) anisotropics, and I will attempt a preliminary comparison with 
theoretical predictions. These notes are intended as a pedagogical introduction 
rather than as a comprehensive review. For comprehehsive discussions of infiation, 
the reader is referred to Refs. 1-3, and for detailed reviews of topological defect 
models to Refs. 4-7. An introduction to quantum field theory methods used in 
modern cosmology can be found in Ref. 8. These notes draw on material presented 
elsewhere^' ^^-^ in which some of the topics are treated in more detail. 

I hope to persuade the reader that cosmology is an exciting area of physics with 
close connections to particle and high energy physics, and with a steady stream of 
new data from astronomy and astrophysics. 

The outline of these lectures is as follows: Section 2 is a review of standard 
cosmology, focusing on its basic principless, its observational support and its prob- 
lems. 

Section 3 is a brief overview of "new cosmology." I argue why, to obtain an 
improved cosmological scenario, we need to treat matter using particle physics 
and field theory. Next, I introduce the basic idea of the infiationary Universe 



scenario and explain how it leads to a solution of some of the problems of standard 
cosmology. In particular, it provides a mechanism for the formation of structure 
in the Universe. The section continues with a brief introduction to the topological 
defect models of structure formation, an explanation for the need of dark matter, 
and a survey of the present models. 

In Section 4, I present the basics of structure formation, beginning with a 
survey of some of the relevant large-scale structure data. Structure in the Universe 
is assumed to grow by gravitational instability. I summarize the essentials of the 
Newtonian theory of cosmological perturbations (valid on length scales smaller than 
the apparent horizon (Hubble radius)) and of the relativistic theory^^-* (required 
to study scales beyond the horizon). I also discuss free streaming. 

Section 5 contains an overview of inflationary Universe models and of the mech- 
anism for the generation and evolution of perturbations which they provide. 

Section 6 presents an overview of topological defect models of structure forma- 
tion. To begin, a classification of defects is given. Next it is shown that in models 
which admit topological defects, they are inevitably formed during a symmetry 
breaking phase transition^^-* . Cosmic string and global texture models of structure 
formation are discussed in detail. In particular, it is pointed out that if defects are 
responsible for seeding galaxies, there must be new physics at a scale of ?] ~ 10^^ 
GeV. 

Section 7 focuses on CMB anisotropies. It is shown why theories of struc- 
ture formation inevitably produce such anisotropies, the predictions of the various 
models are reviewed, and a comparison with recent observations is given. 

In this writeup, units in which c = h = kg = 1 are used unless mentioned oth- 
erwise. The space-time metric Qf^i, is taken to have signature (+, — , — , — ). Greek 
indices run over space and time, latin ones over spatial indices only. The Hub- 
ble expansion rate is H{t) = d{t)/a{t), with a{t) the scale factor of a Friedmann- 
Robert son- Walker (FRW) Universe. The present value of iJ is lOO/i kms^^ Mpc~^, 
where 0.4 < h < 1. Unless stated otherwise, the value of h is taken to be 0.5. The 



cosmological redshift at time t is denoted by z{t). As usual, the symbols G and 
rUpi stand for Newton's constant and Planck mass, respectively. Distances are 
measured in pc ("parsec") or Mpc, where 1 pc corresponds to 3.1 light years. 

2. Review of Standard Cosmology 

2.1. Principles 

The standard big bang cosmology rests on three theoretical pillars: the cosmo- 
logical principle, Einstein's general theory of relativity and a perfect fluid descrip- 
tion of matter. 

The cosmological principle^^-* states that on large distance scales the Universe 
is homogeneous. From an observational point of view this is an extremely non- 
trivial statement. On small scales the Universe looks rather inhomogeneous. The 
inhomogeneities of the solar system are obvious to everyone, and even by the naked 
eye it is apparent that stars are not randomly distributed. They are bound into 
galaxies, dynamical entities whose visible radius is about 10 pc. Telescopic ob- 
servations show that galaxies are not randomly distributed, either. Dense clumps 
of galaxies can be identified as Abell clusters. In turn, Abell cluster positions are 
correlated to produce the large-scale structure dominated by sheets (or filaments), 
with typical scale 100 Mpc, observed in recent redshift surveys^^-*. Until recently, 
every new survey probing the Universe to greater depth revealed new structures 
on the scale of the sample volume. In terms of the visible distribution of matter 
there was no evidence for large-scale homogeneity. This situation changed in 1992 
with the announcement •* that a new redshift survey, complete to a depth of about 
500/i"^ Mpc, had discovered no prominent structures on scales larger than lOOh^^ 
Mpc. This is the first observational evidence from optical measurements in favor 
of the cosmological principle. However, to put this result in perspective we must 
keep in mind that the observed isotropy of the CMB temperature -* to better than 
10~^ on large angular scales has been excellent evidence for the validity of the 
cosmological principle. 



The second theoretical pillar is general relativity, the theory which determines 
the dynamics of the Universe. According to the cosmological principle, space at 
any time t is a three dimensional surface with maximal symmetry (translations 
and rotations). There are three families of such spaces^^-*: flat Euclidean space R^, 
the three sphere S^, and the hypersphere H^. The proper distance ds'^ on these 
three surfaces can be written in spherical coordinates as 



ds"^ = a{tf 



dr^ 



kr' 



+ r2(rf^2^sin2^V 



(2.1) 



The constant k is +1, 0, or —1 respectively for S^, R^ and H^. 

The Einstein equations of general relativity imply that a{t) - called the scale 
factor of the Universe - evolves in time. The proper distance/time in space-time 
is 

dr'^ 



ds^ = dr - a{ty 



+ r^{d^^ + sm^Mip^ 



2.2 



1 — kr"^ 
By a coordinate transformation, a{t) can be set equal to 1 at the present time to- 



Figure 1: Sketch of the expanding Universe. Concentric circles indicate space 
at fixed time, with time increasing as the radius gets larger. Points at rest have 
constant comoving coordinates. Their world lines are straight lines through the 
origin (e.g. L). 



To obtain a simple visualization of an expanding Universe, consider space to 
be the surface of a balloon. We draw a grid on the surface and use it to define 



coordinates xf^' (the superscript c stands for conioving). Points at rest on the surface 
of the balloon have constant comoving coordinates. However, if the balloon is being 
inflated, then the physical distance Ax^ betwen two points at rest with comoving 
separation lS.x^ increases: 

Ax^ = a(t)Aa:^ (2.3) 

The scale factor a(t) is proportional to the radius of the balloon (see Fig. 1). 

According to Einstein's equivalence principle, particles in the absence of exter- 
nal nongravitational forces move on geodesies, curves which extremize ds^ . The 
velocity of a particle relative to the expansion of the Universe is called peculiar 
velocity Vp 

diX 
Vp = ait) — (2.4) 



and obeys the equation 



from which it follows that 



Vp + -vp = 0, (2.5) 



vpit) r^ a-\t) . (2.6) 

The dynamics of an expanding Universe is determined by the Einstein equa- 
tions, which relate the expansion rate to the matter content, specifically to the 
energy density p and pressure p. For a homogeneous and isotropic Universe, they 
reduce to the Friedmann-Robertston- Walker (FRW) equations 

f a\ k SnG , „, 

- -- = ^rP 2.7 

\a ) a^ 3 

a AnG , ^ . ,^ ^. 

- = ^ P + 3p . 2.8 

a 3 

These equations can be combined to yield the continuity equation (with Hubble 



constant H = a/a) 

p = -3H{p + p). (2.9) 

The third key assumption of standard cosmology is that matter is described 
by an ideal gas with an equation of state 

p = wp. (2.10) 

For cold matter, pressure is negligible and hence w = 0. From (2.9) it follows that 

Prnit) ^ a-^t) , (2.11) 

where pm is the energy density in cold matter. For radiation we have w = 1/3 and 
hence it follows from (2.9) that 

Pr{t) r^ a-\t) , (2.12) 

Pr{t) being the energy density in radiation. 
2.2. Observational Pillars 

The first observational pillar of standard cosmology is Hubble's redshift-distance 
relationship^®-' (Fig- 2) 

z = Hd, (2.13) 

where H is the present Hubble expansion constant, d is the distance to a galaxy, 
and z is its redshift 

^-^-1, (2.14) 

Ae(Ao) being the wavelength of light at the time of emission (detection). 



Figure 2: A recent redshift-distance plot of galaxies •* . The distances are deter- 
mined using the Tully-Fisher method. See Ref. 31 for a detailed discussion of the 
method and errors. 

There is an easy intuitive derivation of this result. A wave in an expanding 
background will have a wavelength which increases as the scale factor a{t). Hence 
for light emitted at time te 

.(te) = "-^-l. (2.15) 

a{te) 

For light emitted close to the present time we can Taylor expand the above result 
to obtain 

z{te)c^^Ato-te)^H{to)d. (2.16) 

Equation (2.15) defines the cosmological redshift, which can be used as a measure 
of cosmic time. 

The second observational pillar of standard cosmology is the existence and 
black body nature of the CMB20'21) 



To understand the connection^^) , consider matter in an expanding Universe. 
As we go backwards in time, the density of matter increases as a~^{t), and as a 
consequence the temperature grows. Above a temperature of 13.6 eV, atoms are 
ionized, and a bath of photons in thermal equihbrium must be present. Photons 
still scatter frequently below this temperature. At some time tree the scattering 
length of a photon becomes longer than the Hubble radius. After that, photons 
travel without scattering. At tree, the distribution of photons is of black body type. 
A special feature of black body spectra is that the spectral shape is maintained 
even after tree- The only change is that the temperature redshifts 

nt) = ^^Titree). (2.17) 

Hence, the standard Big Bang model predicts a black body spectrum of photons 
with temperature 

-^o ^ -i-reeZytree) i l^.ioj 

where Tree is the temperature at tree, determined by comparing the largest rate of 
scattering of photons below recombination, that due to Thomson scattering, with 
the Hubble expansion rate, yielding the result 

Tree ^ 0.25 eV~4000°K. (2.19) 

The corresponding redshift is determined by measuring Tq. 

In 1965, Penzias and Wilson^^) discovered this remnant black body radiation at 
a temperature of about 3° K. Since the spectrum peaks in the microwave region it is 
now called the CMB (cosmic microwave background). Recent satellite (COBE)^^-* 
and rocket^"^-* experiments have confirmed the black body nature of the CMB to 
very high accuracy. The temperature is 2.73° K= Tq which corresponds to 

Zitree) = ^ree ~ 10^ . (2.20) 

Given the existence of the CMB, we know that matter has two components: 
dust (with energy density Pm{t)) and radiation (with density Pr{t))- At the present 



time to, Pm{t) ^ Pr{t)- The radiation energy density is determined by Tq, and the 
matter energy density can be estimated by analyzing the dynamics of galaxies 
and clusters and using the virial theorem. However, since by (2.11) and (2.12) 
Pm{t) ~ a{t)~^ and Pr{t) ~ (^{t)~'^i as we go back in time the fraction of energy 
density in radiation increases, and the two components become equal at a time teg, 
the time of equal matter and radiation. The corresponding redshift is 



where 



2eg ~ fi /igo^ io4 (2.21) 



fi = - (to) , (2.22) 

Pc 



Pc being the density for a spatially flat Universe (the critical density), and /150 is 
the value of H in units of 50 km s^^ Mpc^^. 

The time teq is important for structure formation. As we will see in Section 4, it 
is only after tgg that perturbations on scales smaller than the Hubble radius H^^{t) 
can grow. Before then, the radiation pressure prevents growth. A temperature- 
time plot of the early Universe is sketched in Fig. 3, Note that tgg < tree- 



Figure 3: Temperature-time diagram of standard big bang cosmology. The 
present time, time of last scattering and time of equal matter and radiation are 
to, tree and teq respectively. The Universe is radiation-dominated before teq (Re- 
gion A) and matter-dominated in Region B. Before and after treo respectively, the 
Universe was opaque and transparent, respectively, to microwave photons. 
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The third observational pillar of standard big bang cosmology concerns nucleosynthesis^^'^^-*! 
- the production of light elements (heavy elements are formed in supernovae). 
Above a temperature of about 10^ K, the nuclear interactions are sufficiently fast 
to prevent neutrons and protons from fusing. However, below that temperature, 
it is thermodynamically favorable for neutrons and protrons to fuse and form deu- 
terium, helium 3, helium 4 and lithium 7 through a long and interconnected chain 
of reactions. The resulting light element abundances depend sensitively on the 
expansion rate of the Universe and on Qb, the fraction of energy density pB at 
present in baryons relative to the critical density pc- In Fig. 5, recent theoreti- 
cal calculations^^-* of the abundances are shown and compared with observations. 
Demanding agreement with all abundances leaves only a narrow window 

3 X 10~^° < ?7 < 10"^ (2.23) 

where rj is the ratio of baryon number density ub to photon number density n^ 

7] = ^. (2.24) 

From (2.23), it follows that Qb is constrained: 

0.01 <f]B/i2 < 0.035. (2.25) 

In particular, if the Universe is spatially flat and the cosmological constant is 
negligible, there must be nonbaryonic dark matter. We will return to the dark 
matter issue in Section 3. 

The final pillar of standard cosmology is the near isotropy of the CMB^^-*. 
After subtracting the dipole anisotropy which is presumed to be due to the motion 
of the earth relative to the rest frame defined by the CMB, no anisotropies have 
been detected to a level of better than 10~^, i.e., the temperature difference ST{'d) 
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Figure 5: Light element abundances as a function of tlie baryon to entropy ratio 
1] (from Ref. 27). The sohd curves are the predictions of homogeneous big bang 
nucleosynthesis. The observational limits are indicated on the left vetical axis. 
Theory and observations are only consistent for a narrow range of values of rj. 

between two beams pointing in directions in the sky separated by an angle i} (Fig. 
4) satisfies 



T 



< 10"* (2.26) 



12 



Figure 4: Sketch of a CMB anisotropy experiment. Two radio antennas witli 
beam widtli b collect microwave radiation from points in the sky separated by an 
angle T. The difference in beam intensities is measured. 

on all angular scales i). Here T is the average temperature. 

Until recently, the isotropy of the CMB was the only observational support 
for the cosmological principle. Any inhomogeneities of the Universe on length 
scales comparable to the comoving Hubble radius at tree and larger would generate 
temperature anisotropics by a mechanism discussed in detail in Section 7. Hence, 
the near isotropy of the CMB implies that density fluctuations on large scales must 
have been very small in the early Universe. 

To summarize, the observational pillars of standard cosmology are Bubble's 
redshift-distance relation, the existence and black body nature of the CMB, pri- 
mordial nucleosynthesis, and the isotropy of the CMB. Note, in particular, that 
no tests of big bang cosmology say anything about the evolution of the Universe 
before the time of nucleosynthesis. Note, also, that not all astronomers accept the 
above observations as support of the Big Bang model. For a recent criticism see 
Ref. 28 (and Ref. 29 for a reply to the criticism). 
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2.3. Problems 

Standard Big Bang cosmology is faced with several important problems. Only 
one of these, the age problem, is a potential conflict with observations. The others 
which I will focus on here - the homogeneity, flatness and formation of structure 
problems (see e.g. Ref. 30) - are questions which have no answer within the theory 
and are therefore the main motivation for the new cosmological models which will 
be discussed in the rest of these lectures. 

From the FRW equations (2.7) and (2.8) it is easy to calculate the age of 
the Universe, given the expansion law (2.11) which holds throughout most of the 
history of the Universe. For a spatially flat Universe, the age r depends on the 
expansion rate if, i.e. on the constant h which is in the range 0.4 < h <1: 

r~-10V- (2.27) 

Globular cluster ages have been estimated to lie in the range 12 — 18 x 10^ yr. 
Thus, theory and observations are only consistent li h < 0.55 (see e.g. Ref. 31). 
In an open Universe the problem is less severe. Modern cosmological models do 
not add any insight into this problem since they only modify the evolution of the 
Universe at very early times t <^ teq. 

The final three problems mentioned above, the homogeneity, fiatness and for- 
mation of structure problems, provided a lot of the motivation for the development 
of the infiationary Universe scenario"^^-* and will hence be discussed in detail. 

The horizon problem is illustrated in Fig. 6. As is sketched, the co moving 
region ip{trec) over which the CMB is observed to be homogeneous to better than 
one part in 10^ is much larger than the comoving forward light cone if (tree) at tree, 
which is the maximal distance over which microphysical forces could have caused 
the homogeneity: 

to 

^p\tre 



ldta-\t)c^3to(l-(^-^y] (2.28) 
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ifiUec) I dta-\t)-3tl^hlic. (2.29) 







From the above equations it is obvious that ip{trec) ^ ^f{trec)- Hence, stan- 
dard cosmology cannot explain the observed isotropy of the CMB. 



Figure 6: A space-time diagram (physical distance Xp versus time t) illustrating 
the homogeneity problem: the past light cone ip(t) at the time tree oi last scattering 
is much larger than the forward light cone ^/(t) at tree- 
In standard cosmology and in an expanding Universe, fi = 1 is an unstable 
fixed point. This can be seen as follows. For a spatially flat Universe {Q = 1) 

H^ = ^-Pe, (2.30) 

whereas for a nonflat Universe 

H^ + eT' = ^p, (2.31) 



^ = TZ^- (2-32) 



with 

k_ 

The quantity e is proportional to s ' , where s is the entropy density. Hence, in 
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standard cosmology, e is constant. Combining (2.30) and (2.31) gives 

Pc SttG Pc 

Thus, as the temperature decreases, Vt — l increases. In fact, in order to explain the 
present small value oiVt — l ~ 0{1), the initial energy density had to be extremely 
close to critical density. For example, at T = 10^^ GeV, (2.33) implies 

^^^ ~ 10-50 . (2.34) 

Pc 

What is the origin of these fine tuned initial conditions? This is the flatness problem 
of standard cosmology. 

The last problem of the standard cosmological model I will mention is the 
"formation of structure problem." Observations indicate that galaxies and even 
clusters of galaxies have nonrandom correlations on scales larger than 50 Mpc (see 
e.g. Ref. 14). This scale is comparable to the comoving horizon at teg- Thus, if 
the initial density perturbations were produced much before teq-, the correlations 
cannot be explained by a causal mechanism. Gravity alone is, in general, too weak 
to build up correlations on the scale of clusters after teq (see, however, the explosion 
scenario of Ref. 32). Hence, the two questions of what generates the primordial 
density perturbations and what causes the observed correlations, do not have an 
answer in the context of standard cosmology. This problem is illustrated by Fig. 
7. 

Finally, let us address the cosmological constant problem. All known symme- 
tries of nature and principles of general relativity allow for the presence of a term in 
the Einstein equations which acts like matter with energy A and pressure —A, i.e., 
with an equation of state p = —p. If it is not to dominate the present expansion 
rate of the Universe, the cosmological constant A must be very small 

A < 3i7^ ~ lO-^^GeV^ . (2.35) 

On dimensional grounds, we would expect A to be of the order m^, ~ 10"^® GeV . 



'pi 



16 



Figure 7: A sketch (conformal separation vs. time) of the formation of struc- 
ture problem: the comoving separation dc between two clusters is larger than the 
forward light cone at time tgg. 

Thus, the cosmological constant is about 140 orders of magnitude smaller than 
what we would expect it to be (for recent reviews of the cosmological constant 
problem, see Ref. 33). 

As we will see, modern cosmology does not address the cosmological constant 
problem. If anything, the problem will manifest itself in a more apparent manner. 

Due to the formation of structure problem, there can be no causal physical the- 
ory for the origin of structure (with nontrivial spatial correlations) in the Universe 
in the context of the Standard Big Bang theory. The main breakthrough of modern 
cosmology is that it provides solutions to this problem. The key to understand- 
ing this breakthrough in cosmology is the realization of the internal inconsistency 
of the standard picture when extrapolated to times much before nucleosynthesis. 
Standard cosmology is based on the assumption that matter continues to be de- 
scribed by an ideal radiation gas to arbitrarily high temperatures. This is clearly 
in contrast to what nuclear and particle physics tells us. As we go backwards in 
time towards the Big Bang, nuclear physics and eventually particle physics effects 
will take over. To describe matter correctly, a quantum field theoretic description 
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must be used. 

3. New Cosmology and Structure Formation 

The goal of this section is to present an overview of what can be gained if 
we go beyond standard cosmology and allow matter to be described in terms of 
concepts from particle physics. Detailed discussions of the models will be given in 
later sections. 

3.1. The Inflationary Unvierse 

The idea of inflation^^-* is very simple. We assume there is a time interval 
beginning at ti and ending at t^ (the "reheating time" ) during which the Universe 
is exponentially expanding, i.e., 

a(t)~e■^^ te[ti,tR] (3.1) 

with constant Hubble expansion parameter H. Such a period is called "de Sitter" 
or "inflationary." The success of Big Bang nucleosynthesis sets an upper limit to 
the time of reheating: 

tR < tNS , (3.2) 

tNS being the time of nucleosynthesis. 



Figure 8: The phases of an inflationary Universe. The times U and Ir denote 
the beginning and end of inflation, respectively. In some models of inflation, there 
is no initial radiation domintated FRW period. Rather, the classical space-time 
emerges directly in an inflationary state from some initial quantum gravity state. 
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The phases of an inflationary Universe are sketched in Fig. 8. Before the onset 
of inflation there are no constraints on the state of the Universe. In some models a 
classical space-time emerges immediately in an inflationary state, in others there is 
an initial radiation dominated FRW period. Our sketch applies to the second case. 
After t^, the Universe is very hot and dense, and the subsequent evolution is as 
in standard cosmology. During the inflationary phase, the number density of any 
particles initially in thermal equilibrium aX t = ti decays exponentially. Hence, the 
matter temperature Tm(t) also decays exponentially. At t = tR, all of the energy 
which is responsible for inflation (see later) is released as thermal energy. This is 
a nonadiabatic process during which the entropy increases by a large factor. The 
temperature-time evolution in an inflationary Universe is depicted in Fig. 9. 



Figure 9: The time dependence of matter temperature in an inflationary Universe. 
During the period of exponential expansion, the temperature decreases exponen- 
tially. At the end of inflation the energy density of the scalar fleld responsible 
for inflation is transferred to ordinary matter. This leads to reheating. The criti- 
cal temperature Tc is the temperature at which the initial matter thermal energy 
density becomes less than the scalar fleld energy density (see Chapter 5). 

Fig. 10 is a sketch of how a period of inflation can solve the homogeneity 
problem. At = t^i — ti is the period of inflation. During inflation, the forward 
light cone increases exponentially compared to a model without inflation, whereas 
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the past light cone is not affected for t > Ir. Hence, provided At is sufficiently 
large, if{tR) will be greater than ipit^). The condition on At depends on the tem- 
perature Tji corresponding to time t^, the temperature of reheating. Demanding 
that ifitR) > ipitRJ we find, using the analogs of (2.28) and (2.29), the following 
criterion 



for Tji ~ 10^'^GeV and Tq ~ 10~^"^GeV (the present microwave background tem- 
perature). Thus, in order to solve the homogeneity problem, a period of inflation 
with 

At>50iJ"^ (3.4) 

is required. 

Inflation also can solve the flatness problem^^'"^^-* The key point is that the 
entropy density s is no longer constant. As will be explained later, the temperatures 
at ti and t^ are essentially equal. Hence, the entropy increases during inflation 
by a factor exp(3i7At). Thus, e decreases by a factor of exp(— 2/JAt). With the 
numbers used in (3.3): 

Cafter ~ 10 Cbefore • (3-5) 

Hence, (p — Pc)/p can be of order 1 both at ti and at the present time. In fact, if 
inflation occurs at all, the theory then predicts that at the present time fi = 1 to 
a high accuracy (now ^2 < 1 would require special initial conditions). 

What was said above can be rephrased geometrically: during inflation, the 
curvature radius of the Universe - measured on a fixed physical scale - increases 
exponentially. Thus, a piece of space looks essentially flat after inflation even if it 
had measurable curvature before. 

Most importantly, inflation provides a mechanism which in a casual way gen- 
erates the primordial perturbations required for galaxies, clusters and even larger 
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Figure 10: Sketch (physical coordinates vs. time) of the solution of the homo- 
geneity problem. During inflation, the forward light cone lf{t) is expanded expo- 
nentially when measured in physical coordinates. Hence, it does not require many 
e-foldings of inflation in order that lf{t) becomes larger than the past light cone 
at the time of last scattering. The dashed line is the forward light cone without 
inflation. 

objects. In inflationary Universe models, the Hubble radius ("apparent" horizon), 
3t, and the "actual" horizon (the forward light cone) do not coincide at late times. 
Provided (3.3) is satisfied, then (as sketched in Fig. 11) all scales within our appar- 
ent horizon were inside the actual horizon since ti. Thus, it is in principle possible 
to have a casual generation mechanism for perturbations^^^^^-* . 

The generation of perturbations is supposed to be due to a causal microphysical 
process. Such processes can only act coherently on length scales smaller than the 
Hubble radius init) where 

iH{t)=H-\t). (3.6) 
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Figure 11: A sketch (physical coordinates vs. time of the solution of the formation 
of structure problem. Provided that the period of inflation is sufficiently long, the 
separation dc between two galaxy clusters is at all times smaller than the forward 
light cone. The dashed line indicates the Hubble radius. Note that dc starts out 
smaller than the Hubble radius, crosses it during the de Sitter period, and then 
reenters it at late times. 

A heuristic way to understand the meaning of ^/^(t) is to realize that it is the 
distance which light (and hence the maximal distance any causal effects) can prop- 
agate in one expansion time: 

£jy(t)~a(t) I a{t')-^dt'. (3.7) 

t 

In Section 5 a more mathematical justification for the definition and role of ^_H"(t) 
will be given. 
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As will be discussed in Section 5, the density perturbations produced during in- 
flation are due to quantum fluctuations in the matter and gravitational flelds"^^'^^-* . 
The amplitude of these inhomogeneities corresponds to a tempertuare Th 

Th-^H, (3.8) 

the Hawking temperature of the de Sitter phase. This implies that at all times 
t during inflation, perturbations with a flxed physical wavelength ~ H^^ will be 
produced. Subsequently, the length of the waves is streched with the expansion 
of space, and soon becomes larger than the Hubble radius. The phases of the 
inhomogeneities are random. Thus, the inflationary Universe scenario predicts 
perturbations on all scales ranging from the comoving Hubble radius at the be- 
ginning of inflation to the corresponding quantity at the time of reheating. In 
particular, provided that inflation lasts sufficiently long (see (3.4)), perturbations 
on scales of galaxies and beyond will be generated. 

Now that the reader is (hopefully) convinced that inflation is a beautiful idea, 
the question arises how to realize this scenario. The initial hope was that the 
same scalar flelds (Higgs flelds) which particle physicists introduce in order to 
spontaneously break the internal symmetries of their fleld theory models would 
lead to inflation. This hope was based on the fact that the energy density p and 
pressure p of a real scalar fleld (^(x, t) are given by 



p{v) = \^^ + \{V^f + V{^) 



(3.9) 



Thus, provided that at some initial time tj 

^{x,ti)=V^{xiti) = Q (3.10) 

and 

nv^feti))>0, (3.11) 
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the equation of state of matter will read 

P=-P (3.12) 

and, from the FRW equations it will follow that 



a(t) = e^^ , if2 = ^\/(v9). (3.13) 



The next question is how to realize the required initial conditions (3.10) and 
to maintain the key constraints 

V92 « \/(v9) , (Vvp)2 « \/(v9) (3.14) 

for sufficiently long (see (3.4)). This typically requires potentials which are very 
flat near (/^(a;, tj). Worse yet, the parameters of the potential V{ip) must be chosen 
such that the final amplitude of density perturbations is sufficiently small to agree 
with the limits on CMB anisotropics. As we will demonstrate in Section 5, these 
conditions impose severe constraints on the constants which appear in V{ip). 

In light of these difficulties it is important to keep in mind that inflation can 
also be generated by modifying gravity at high curvatures (see e.g., Refs. 39-41). 
It is also wise to investigate alternative theories of structure formation which do 
not rely on inflation. 

3.2. Topological Defect Models 

According to particle physics theories, matter at high energies and tempera- 
tures must be described in terms of flelds. Gauge symmetries have proved to be 
extremely useful in describing the standard model of particle physics, according to 
which at high energies the laws of nature are invariant under a nonabelian group 
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G of internal symmetry transformations 

G = SU(3)c X SU(2)l X U{1)y (3.15) 

which at a temperature of about 200 MeV is spontaneously broken down to 

G" = SU(3)c X U(l) . (3.16) 

The subscript on the SU(3) subgroup indicates that it is the color symmetry group 
of the strong interactions, SU(2)7;,x U(l)y is the Glashow-Weinberg-Salam (WS) 
model of weak and electromagnetic interactions, the subscripts L and Y denoting 
left handedness and hypercharge respectively. At low energies, the WS model 
spontaneously breaks to the U(l) subgroup of electromagnetism. 

Spontaneous symmetry breaking is induced by an order parameter ip taking 
on a nontrivial expectation value < ip > below a certain temperture T^.. In some 
particle physics models, v? is a fundamental scalar field in a nontrivial representa- 
tion of the gauge group G which is broken. However, ip could also be a fermion 
condensate. 

The transition taking place at T = Tg is a phase transition and T^ is called 
the critical temperature. From condensed matter physics it is well known that in 
many cases topological defects form during phase transitions, particularly if the 
transition rate is fast on a scale compared to the system size. When cooling a metal, 
defects in the crystal configuration will be frozen in; during a temperature quench 
of ^He, thin vortex tubes of the normal phase are trapped in the superfluid; and 
analogously in a temperature quench of a superconductor, flux lines are trapped 
in a surrounding sea of the superconducting Meissner phase (see Fig. 12 for an 
example^^-* of defect formation). 

In cosmology, the rate at which the phase transition proceeds is given by 
the expansion rate of the Universe. Hence, topological defects will inevitably be 
produced in a cosmological phase transition^^-' , provided the underlying particle 
physics model allows such defects. 
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Figure 12: A simulation of defect formation in the 2 + 1 dimensional Abelian 
Higgs model**^-* in an expanding background. The total energy density is plotted 
against the two spatial coordinates. The initial conditions were specified by thermal 
initial conditions with random phases of the order parameter on Hubble scales. At 
later times, the thermal noise has redshifted away, leaving behind trapped energy 
density in vortices. 

Topological defects can be point-like (monopoles), string-like (cosmic strings) •* 
or planar (domain walls), depending on the particle physics model (see Section 6). 
Also of importance are textures^^'^^-* , point defects in space-time. 

Topological defects represent regions in space with trapped energy density. 
These regions of surplus energy can act as seeds for structure formation as is 
illustrated in Fig. 13. For point-like defects, the force which causes clustering 
about the seed can be understood using Newtonian gravity. The process is called 
gravitational accretion. For precise calculations (and in the case of other defects), 
general relativistic effects must be taken into account (see Section 6). 
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Figure 13: Sketch of the basic gravitational accretion mechanism. The topological 
defect (in this case a cosmic string loop) is a configuration of trapped energy- 
density. This excess density produces a Newtonian gravitational attractive force 
on the surrounding matter. 

No stable topological defects arise in the breaking of the WS model. However, 
there is good evidence for phase transitions at very high energies. The coupling 
constants of SU(3)c, SU(2) and U(l) are seen to converge at an energy scale rj of 
about 

T^-lO^^GeV. (3.17) 

It is therefore not unreasonable to speculate that the standard model results from 
the breaking of a larger symmetry group Go at a scale rj. A large class of unified 
gauge theories"^^-* based on a symmetry breaking 

Go — ^G = SU(3) X SU(2) xU(l) (3.18) 

admit topological defects, many theories have cosmic string solutions. 
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Topological defect models of structure formation will be discussed in detail 
in Section 6. Here I will briefly point out by which mechanism correlations on 
all cosmological scales are induced. To be concrete, I consider the cosmic string 
model^^'^®-* . Cosmic strings are one-dimensional defects without ends. Hence, they 
must be either infinite in length or else closed loops. The fact that at the time of 
the phase transition tc, the order parameter has random phases on scales larger 
than the initial correlation length implies that at tc a random walk-like network 
of infinite strings will form^^-*. This implies nontrivial correlations of structures 
seeded by these strings on all scales larger than the initial correlation length. 

3.3. Need for Dark Matter 

At this point we have illustrated two classes of mechanisms by which structure 
formation in the Universe can be seeded: quantum fluctuations during a period 
of inflation, and topological defects. To completely specify a theory of structure 
formation, however, we must also specify what the "dark matter" which dominates 
the energy density of the Universe today consists of (see e.g. Ref. 49). 

The evidence for dark matter has been accumulating over the past decade. 
Measurements of galaxy velocity rotation curves indicate that a large fraction of 
the mass of a galaxy does not shine. In Fig. 14, the measured rotation velocity v is 
plotted as a function of the distance r from the center of the galaxy. This velocity 
is compared to the velocity v{r) which results from the virial theorem, assuming 
that light traces mass (the luminosity profile is given in the upper frame). The 
data is for the spiral galaxies NGC2403 and NGC3198^'^'*. The comparison shows 
that the mass of these galaxies extends significantly beyond the visible radius and 
that - assuming Newtonian gravity is applicable - a large fraction of the mass of 
a spiral galaxy must be dark. 

The fraction of matter that shines can be expressed as a fraction Q\um of the 
critical density pc- The present estimates give 

^lum<0.01, (3.19) 
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Figure 14: Velocity rotation curves for two galaxies (from Ref. 50). The upper 
panel shows the luminosity of the galaxy as a function of the distance from the 
center, the lower panel presents the velocity (vertical axis, in units of kms^ ) as 
a function of the same distance. The solid curve is the velocity inferred from 
the observed luminosity curve, using the virial theorem, the dotted curves are the 
observational results. The fact that the velocity rotation curves remain constant 
beyond the visible radius of the galaxy is strong evidence for galactic dark matter. 

whereas the fraction Vtg of mass in galaxies is much larger^^-* 

f]g~0.03. (3.20) 

It is also possible to estimate the fraction Vtcl of mass which is gravitationally 
bound in clusters. Current estimates using the virial theorem give ' 

fie/ ~ 0.1 -0.2 (3.21) 
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(this comes mainly from studying the infall of galaxies towards the Virgo cluster). 

Finally, the amount ^lss of mass in large-scale structures can be estimated 
by measuring the large-scale peculiar velocities and inferring the mass required to 
generate such velocities. Current estimates give^^'^"^-* 

^^55 = 0.8 ±0.5. (3.22) 

As mentioned in Section 2, nucleosynthesis provides independent limits on the 
fraction Qg of mass in baryons (see (2.25)). Comparing (3.19-3.22) with (2.25) we 
conclude: 

i) There must be baryonic dark matter. In fact, most of the dark matter in 
galaxies and/or clusters could be baryonic, and some of this baryonic dark 
matter may have recently been discovered by gravitational microlensing. 

ii) liQ = 1, then most of the matter in the Universe consists of nonbaryonic dark 
matter. In fact, there is increasing evidence (see (3.22)) that nonbaryonic 
dark matter must exist independent of the theoretical prejudice for fi = 1. 

The dark matter in the Universe is visible only through its gravitational effects. 
Hence, nonbaryonic dark matter candidates can be divided into two classes, cold 
dark matter (CDM) and hot dark matter (HDM). 

CDM particles are cold, i.e., their peculiar velocity v is negligible at the time 
teq when structure formation begins: 

v{teq) < 1 . (3.23) 

Candidates for CDM include the axion (coherent oscillations of a low mass scalar 
field) and neutralinos (the lightest stable supersymmetric particle, which must be 
neutral). 
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HDM particles are relativistic at teq: 

viteg) ~ 1 . (3.24) 

The prime candidate is a SS/igQ eV tau neutrino. Note that this mass is well within 
the experimental bounds for the tau neutrino mass, and also that many particle 
physics models - in particular those which lead to neutrino oscillations - predict 
masses of this order of magnitude. 

3.4. Survey of Models 

Any theory of structure formation must specify both the source of fluctuations 
and the composition of the dark matter. The reader is warned that the model 
called the "CDM Model" is a model with CDM AND perturbations generated by 
quantum fluctuations during a hypothetical period of inflation. 

Inflation-based models were the first to be considered in quantitative detail, 
initially assuming a HDM-dominated Universe. Almost immediately, however, 
contradictions with basic observations appeared^"^-* . 

The problem of HDM-based inflationary models is related to neutrino free 
streaming^^-* . The primordial perturbations in this theory are dark matter fluc- 
tuations, but because of the large velocity of the dark matter particles, the inho- 
mogeneities are washed out on all scales below the neutrino free streaming length 
A,^(t), 

\]{t) r^v{t)z{t)t, (3.25) 

which is the comoving distance the particles move in one Hubble expansion time. 
Since the neutrino velocity v{t) and the redshift z{t) both scale as a(t)^^, the free 
streaming length decreases as 

after teq (before teq the radiation pressure dominates). Hence, in an inflationary 
HDM model all perturbations on scales A smaller than the maximal value of A^(t) 
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are erased. The critical scale A^^^ is given by the value of A^(t) at the time when 
the neutrinos become non-relativistic which is in turn determined by the neutrino 
mass m,y. The result is 



A— ^SOMpc^^j , (3.27) 

a scale much larger than the mean separation of galaxies and clusters. Since we 
observe galaxies outside of large-scale structures, this model is in blatant disagree- 
ment with observations. 

Inflation-based models are hence only viable if (at least a substantial fraction 
of) the dark matter is cold. Such models have become known as "CDM mod- 
els" , and are to a first approximation rather successful at predicting the clustering 
properties of galaxies and galaxy clusters^^-* . There are many parameters in CDM 
models: the amplitude of the density perturbations, the power of the spectrum 
(see Section 4), the value of fi, the fraction VLb of baryons, to mention some of 
the main ones. It is also possible to add a small fraction Vty of hot dark matter 
(yielding a class of so-called "Mixed Dark Matter" models). 

Topological defect models were first developed in the context of CDM. Theo- 
ries based on cosmic strings^^^^^-* or on global textures^^'^^-* have also been fairly 
successful in explaining observations (again to a first approximation). 

It is important to note that if perturbations are seeded by long-lived topological 
defects (e.g., cosmic strings), then the above arguments against hot dark matter 
disappear^^'^^-*. The seed perturbations can survive neutrino free streaming as long 
as the seeds remain present for many Hubble expansion times. If we consider a 
comoving scale A much smaller than \^^^ of (3.27), then a dark matter perturbation 
will begin to grow about the seed fiuctuations at a time t(A) when 

A^(t(A)) = A. (3.28) 

Cosmic string based hot dark matter models have also been successful at explaining 
the qualitative features of observations^^-* . 
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4. Basics of Structure Formation 

In the structure formation models mentioned in the previous section, small 
amplitude seed perturbations are predicted to arise due to particle physics effects 
in the very early Universe. They then grow by gravitational instability to produce 
the cosmological structures we observe today. In order to be able to make the 
connection between particle physics and observations, it is important to understand 
the gravitational evolution of fluctuations. This section will introduce the basic 
concepts of this topic. We begin, however, with an overview of some of the relevant 
data. 

4.1. Survey of Data 

It is length scales corresponding to galaxies and larger which are of great- 
est interest in cosmology when attempting to find an imprint of the primordial 
fluctuations produced by particle physics. On these scales, gravitational effects 
are assumed to be dominant, and the fluctuations are not too far from the lin- 
ear regime. On smaller scales, nonlinear gravitational and hydrodynamical effects 
determine the final state and mask the initial perturbations. 

To set the scales, consider the mean separation of galaxies, which is about 5h~^ 
Mpc -*, and that of Abell clusters which is around 25h^ Mpc •*. The largest 
coherent structures seen in current redshift surveys have a length of about 100h~^ 
Mpc, the recent detections of CMB anisotropics probe the density field on length 
scales of about 10"^ /i~^ Mpc, and the present horizon corresponds to a distance of 
about 3 • lO^h^^ Mpc. 

Galaxies are gravitationally bound systems containing billions of stars. They 
are non-randomly distributed in space. A quantitative measure of this non-randomness 
is the "two-point correlation function" ^2(^) which gives the excess probability of 
finding a galaxy at a distance r from a given galaxy: 

fe(,) =< !!(n^ > . (4.1) 

no 
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Here, hq is the average number density of galaxies, and n{r) is the density of 
galaxies a distance r from a given one. The pointed braces stand for ensemble 
averaging. 



Figure 15: Recent observational results for the two point correlation function of 
IRAS galaxies, in both a volume limited subsample and a complete galaxy sample 
(see Ref. 66 which explains the meaning of the measure ^3(^)). 

Recent observational results from a redshift survey of IRAS (infrared) galaxies 
yields reasonable agreement^^-* with a form (see Fig. 15) 

6(r) ^ [^y (4.2) 
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with scaling length ro — 5/i Mpc and power 7 ~ 1.8. A theory of structure for- 
mation must explain both the amplitude and the slope of this correlation function. 

Galaxies do not all have the same mass. There are more smaller galaxies than 
large ones (our galaxy is a large one). The distribution of galaxy masses is given 
by the "galaxy mass function" n{M), where n{M)dM is the number density of 
galaxies in the mass range [M, M + dM]. Since we can only measure luminos- 
ity but not mass, the observable measure of the galaxy distribution is (piL), the 
"galaxy luminosity function." Now, (f){L)dL is the number density of galaxies with 
luminosity L in the interval [L, L + dL]. Recent results for (f){L) from the IRAS 
galaxy survey are reproduced in Fig. 16^^-*. 



Figure 16: The luminosity function of IRAS galaxies. The vertical axis is <I){L), 
the horizontal axis is luminosity L in units of solar luminosity. The solid curve 
represents the data, the dashed curve is a fit to a theoretical model (see Ref. 67). 
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Theories must also be able to explain the internal mass distribution of galaxies 
which can be inferred from the galaxy rotation curves (see Fig. 14). According to 
the virial theorem, the velocity v{r) at a distance r from the center of a galaxy is 
determined by the mass M{r) inside of r: 

m.2 ^rnAm (43) 






where m is a test mass. Hence 



„M=,G^V'\ (44) 



A constant rotation velocity implies M{r) ~ r and hence a density profile p{r) of 

p{r) r^r^'^ . (4.5) 

This condition puts constraints on the possible composition of the galactic dark 
matter. 

An Abell cluster^®) is a region in space with greater than fifty bright galaxies 
in a sphere of radius 1.5h^^ Mpc, i.e., a region with a very high overdensity of 
galaxies. Observations indicate that Abell clusters are not distributed randomly 
in space. The cluster two point correlation function ^c(^) has a form similar to 

(4 2)65,69). 

e,(r)^(^)' (4.6) 

with r-Q ~ I5h^^ Mpc and 7 ~ 2 (see Fig. 17 which is taken from a recent analysis 
of rich clusters of galaxies selected from the APM Galaxy Survey^^-*). 

There are several remarkable features about the clustering of galaxies and 
galaxy clusters. First, there is evidence for universality of the functional form of 
^{r); its slope is about 2 for both populations. Secondly, relative to their respective 
mean separations, galaxies are more clustered than galaxy clusters. This can be 
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Figure 17: The two point correlation function C,c (denoted by ^cc in the figure) 
of clusters of galaxies drawn from the APM galaxy survey as a function of their 
separation s. Results for two different samples of clusters are given. 

explained by the action of gravity. Gravity has had longer to act on the scales of 
galaxies than on that of clusters, and has hence amplified the galaxy correlation 
function relative to that of clusters. 

There is a wide spread of cluster masses which can be described by the clus- 
ter multiplicity function $c(^)) where ^(.{n)dn is the number density of clusters 
containing between n and n + dn galaxies. Fig. 18 is a sketch of the observed 
cluster multiplicity function taken from Ref. 70. Note that the cluster mass func- 
tion inferred from Fig. 18 and the galaxy mass function deduced from the galaxy 
luminosity function of Fig. 16 match up quite well at a mass of about IO^^Mq. 
Below this mass, the objects are well defined dynamical entities whereas for larger 
masses they are composed of fragments. A reason for the difference may be due to 
the fact that clouds of more than lO^^M© cannot cool without fragmenting ^ 

On scales larger than galaxy clusters there is not at present a clear mathemati- 
cal description of structure. Many galaxy redshift surveys have discovered coherent 
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Figure 18: The cluster multiplicity function $c(^) (denoted as n on the figure) 
as a function of the number n (horizontal axis). 

filamentary and planar structures and voids on scales of up to lOO/i^^ Mpc^"*'^^^^^-* . 
For example, the astronomers working on the "Center for Astrophysics" redshift 
survey^'*-' have analyzed many adjacent slices of the northern celestial sphere. For 
all galaxies above a limiting magnitude of 15.5 they measured the redshifts z. Fig. 
19 is a sketch of redshift versus angle a in the sky for one slice. The second direction 
in the sky has been projected onto the a — z plane. The most prominent feature is 
the band of galaxies at a distance of about lOO/i^^ Mpc. This band also appears 
in neighboring slices and is therefore presumably part of a planar density enhance- 
ment of comoving planar size of at least (50 x 100) x h^'^ Mpc^. This structure is 
often called the "great wall." It is a challenge for theories of structure formation 
to explain both the observed scale and topology of the galaxy distribution. 

Until 1992 there was little evidence for any convergence of the galaxy distribu- 
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Figure 19: Results from the CFA redshift survey. Radial distance gives the 
redshift of galaxies, the angular distance corresponds to right ascension. The results 
from several slices of the sky (at different declinations) have been projected into 
the same cone. 

tion towards homogeneity. Each new survey led to the discovery of new coherent 
structures in the Universe on a scale comparable to that of the survey. In 1992, 
preliminary results of a much deeper redshift survey were announced^^-* which for 
the first time found no new coherent structures on scales larger than lOOh^^ Mpc. 
This is the first direct evidence for the cosmological principle from optical surveys 
(the isotropy of the CMB has for a long time been a strong point in its support). 

In summary, a lot of data from optical and infrared galaxies alone are currently 
available, and new data are being collected at a rapid rate. The observational 
constraints on theories of structure formation are becoming tighter. A lot of the- 
oretical work is needed in order to allow for detailed comparisons between theory 
and observations. 
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4.2. Gravitational Instability 

In this article we only discuss theories in which structures grow by gravitational 
accretion. The basic mechanism is illustrated in Fig. 20. Consider first a flat 
space-time background. A density perturbation with 5p > will then give rise 
to an excess gravitational attractive force F acting on the surrounding matter. 
This force is proportional to 5p, and will hence lead to exponential growth of the 
perturbation since 

Sp ^ F r^ Sp ^ Sp ^ exp(«t) (4-7) 

with some constant a. 



Figure 20: Sketch of the gravitational instability mechanism. The vertical axis 
is the density perturbation Z^ as a function of a line in space (x). A small initial 
overdensity (A) will cause a gravitational acceleration g towards it, which will lead 
to an increase in the perturbation (B). Note that in general underdense regions 
develop in addition to the growing overdense areas. 

In an expanding background space-time, the acceleration is damped by the 
expansion. If r{t) is the physcial distance of a test particle from the perturbation, 
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then on a scale r 

5p~F~-|^, (4.8) 

which results in power-law increase of 5 p. The goal of this subsection is to discuss 
the growth rates of inhomogeneities in more detail (see e.g. Refs. 76 and 77 for 
modern reviews). 

Because of our assumption that all perturbations start out with a small ampli- 
tude, we can linearize the equations for gravitational fluctuations. The analysis is 
then greatly simplifled by going to Fourier space in which all modes 5{k) decouple. 
We expand the fractional density contrast 5{x) as follows: 

^ M^ ^ (27r)-3/VV2 / d'k e'^-^5{k) , (4.9) 

P J 

where K is a cutoff volume which disappears from all physical observables. 
The "power spectrum" P{k) is deflned by 

P{k)=<\5{k)\^>, (4.10) 

where the braces denote an ensemble average (in most structure formation models, 
the generation of perturbations is a stochastic process, and hence observables can 
only be calculated by averaging over the ensemble. For observations, the braces 
can be viewed as an angular average). 

The physical measure of mass fluctuations on a length scale A is the r.m.s. mass 
fluctuation 6M/M{X) on this scale. It is determined by the power spectrum in the 
following way. We pick a center Xq of a sphere B\{xo) of radius A and calculate 

1^1 (a. ^) = l J ^^^(^)^7(^l ' (4.11) 

Bx{xo) 

where V{B\) is the volume of the sphere. Inserting the Fourier decomposition 
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(4.3) we obtain 

,(5M,2. .. V 1 



a, A) =-^ — i-2 j d^xi j d''x2 j d%d%e'^^^-^^y^' 



M 

Bx{0) Bx{0) 



e»fci-£ie-»^2-£25(^^)*5(^ 



2) 



(4.12) 



Taking the average value of this quantity over all Xq yields 



6M^ ^ 



< ( ^ ) (A) >= / d'kWkiXMkr (4.13) 



with a window function Wk{X) with the following properties 



, - 1 k<kx = 1txI\ 
W^{\){ ^ , / (4.14) 

i k> kx ■ 



Therefore the r.m.s. mass perturbation on a scale A becomes 



<\-J^W\ >-klP{kx). (4.15) 



Astronomers usually assume that P{k) grows as a power of k: 

P{k)r^k'', (4.16) 

n being called the index of the power spectrum. For n = 1 we get the so-called 
Harrison-Zel'dovich scale invariant spectrum^^-* . 

Both inflationary Universe and topological defect models of structure formation 
predict a roughly scale invariant spectrum. The distinguishing feature of this 
spectrum is that the r.m.s. mass perturbations are independent of the scale k 
when measured at the time ti{{k) when the associated wavelength is equal to the 
Hubble radius, i.e., when the scale "enters" the Hubble radius. Let us derive this 
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fact for the scales entering during the matter dominated epoch. The time t//(/c) is 
determined by 

k-'a{tH{k)) = tH{k) (4.17) 

which leads to 

tH{k)^k^^. (4.18) 

According to the linear theory of cosmological perturbations discussed in the fol- 
lowing subsection, the mass fluctuations increase as a{t) for t > teg- Hence 

— {k, tnm = y-^j -Jf (k, t) ~ const , (4.19) 

since the first factor scales as k~'^ and - using (4.15) and inserting n = 1 - the 
second as /c^. 

4.3. Newtonian Theory of Cosmological Perturbations 

The Newtonian theory of cosmological perturbations is an approximate analysis 
which is valid on wavelengths A much smaller than the Hubble radius t and for 
negligible pressure p, i.e., p <^ p. It is based on expanding the hydrodynamical 
equations about a homogeneous background solution. 

The starting points are the continuity, Euler and Poisson equations 

P + Y.ipv) = (4.20) 

v + {v-'V)v = -'V(f)--'Vj) (4.21) 

VV = AttGp (4.22) 

for a fluid with energy density p, pressure p, velocity v_ and Newtonian gravitational 
potential 0, written in terms of physical coordinates (t, r). 
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The transition to an expanding space is made by introducing conioving coor- 
dinates x_ and peculiar velocity m = x: 

r = a{t)x (4.23) 

v = a(t)x + a{t)u. (4.24) 

The first term on the right hand side of (4.24) is the expansion velocity. 

The perturbation equations are obtained by linearizing Equations (4.20-22) 
about a homogeneous background solution p = p{t), p = and u = 0. Using the 
definition 

i^'-^, (4.25) 

P 

the linearization ansatz can be written 

p{x,t)=p{t){l + 6{x,t)). (4.26) 

If we consider adiabatic perturbations (no entropy density variations), then after 
some algebra the linearized equations become 

^ + V ■ M = , (4.27) 

u + 2-u = -a^(VS(f) + ciVS) (4.28) 

a 

and 

V^S(j) = AnGpa^S , (4.29) 

with the speed of sound Cg given by 

c5 ^ I . ,4.30) 

The two first order equations (4.27) and (4.28) can be combined to yield a single 
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second order differential equation for 6. Witli tlie lielp of (4.29) tliis equation reads 



,2 



c: 



6 + 2H6-47rGp6-^V^6 = (4.31) 



a 



wliicli in Fourier space becomes 



c2^2 
2 



Sk + 2H6k+ {^-4:7rGp]Sk = 0. (4.32) 



Here, H{t) as usual denotes the expansion rate, and 6k stands for 6{k). 

Already a quick look at Equation (4.32) reveals the presence of a distinguished 
scale for cosmological perturbations, the Jeans length 

A, = ^ (4.33) 

with 

On length scales larger than Aj, the spatial gradient term is negligible, and the 
term linear in 6 in (4.32) acts like a negative mass square quadratic potential with 
damping due to the expansion of the Universe, in agreement with the intuitive 
analysis leading to (4.7) and (4.8). On length scales smaller than Aj, however, 
(4.32) becomes a damped harmonic oscillator equation and perturbations on these 
scales decay. 

For t > teq and for A ^ Aj, Equation (4.32) becomes 

h + ^t^k-^Sk = (4.35) 

and has the general solution 

4(t) = cit2/3 + c2t-i. (4.36) 

This demonstrates that for t > teq and A ^ Aj, the dominant mode of pertur- 
bations increases as a{t), a result we already used in the previous subsection (see 
(4.19)). 
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For A <^ Aj and t > teq, Equation (4.32) becomes 

Sk + 2H6k + c'i(^\ 6k = 0, (4.37) 

and has solutions corresponding to damped oscillations: 

6k{t) ~ a-^/^{t) exp{±icsk / dt'a{t'y'^} . (4.38) 

As an important application of the Newtonian theory of cosmological perturba- 
tions, let us compare sub-horizon scale fluctuations in a baryon-dominated Universe 
{Q = Qb = 1) and in a CDM-dominated Universe with ^c'DM = 0.9 and Q = 1. 
We consider scales which enter the Hubble radius at about teq- 

In the initial time interval teq < t < tree, the baryons are coupled to the 
photons. Hence, the baryonic fluid has a large pressure pB 



PB-Pr = -pr- (4.39) 

Hence, the speed of sound is relativistic 

I!L.y'\j_(£L.V'\ (4.40) 

PmJ V3 \Pm/ 

The value of Cs slowly decreases in this time interval, attaining a value of about 
1/10 at tree- From (4.34) it follows that the Jeans mass Mj, the mass inside a 
sphere of radius Aj, increases until tree when it reaches its maximal value MJ*"^ 

Mr" = Mj{trec) = ^ \Atrecfp{trec) ~ lO^\nhY^/^MQ . (4.41) 



At the time of recombination, the baryons decouple from the radiation fluid. 
Hence, the baryon pressure pB drops abruptly, as does the Jeans length (see (4.34)). 
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The remaining pressure pB is determined by the temperature and thus continues to 
decrease as t increases. It can be shown that the Jeans mass continues to decrease 
after tree, starting from a value 

MjiUec) --lO'^i^h^y^/^MQ (4.42) 

(where the superscript "— " indicates the mass immediately after teq- 

In contrast, CDM has negligible pressure throughout the period t > tgq and 
hence experiences no Jeans damping. A CDM perturbation which enters the Hub- 
ble radius at teq with amplitude 6i has an amplitude at tree given by 

Ok V-rec) ^ -TTT ^^ ' ^^■^^> 

'J'V^eq) 

whereas a perturbation with the same scale and initial amplitude in a baryon- 
dominated Universe is damped 

In order for the perturbations to have the same amplitude today, the initial size of 
the inhomogeneity must be much larger in a BDM-dominated Universe than in a 
CDM-dominated one: 

^r^M-f#4V^'^.^''''M- (4.45) 



\^\teqj / 

For Q = 1 and h = 1/2 the enhancement factor is about 30. 

In a CDM-dominated Universe the baryons experience Jeans damping, but 
after tree the baryons quickly fall into the potential wells created by the CDM 
perturbations, and hence the baryon perturbations are proportional to the CDM 
inhomogeneities (see Fig. 21). 

47 



Figure 21: Comparison of tlie growtli of fluctuations in baryon and CDM dom- 
inated Universes. The horizontal axis is time, the vertical axis is the density 
perturbation. The curve labelled by a describes the evolution of fluctuations in a 
BDM Universe. The growth of CDM perturbations follows the curve 6, and curve c 
is a sketch of the time development of baryon perturbations in a CDM dominated 
Universe. To be specific, perturbations on a scale which enters the Hubble radius 
at teq are considered. 

The above considerations, coupled with information about CMB anisotropies, 
can be used to rule out a model with Vt = Qb = 1- The argument goes as follows 
(see Section 7). For adiabatic fluctuations, the amplitude of CMB anisotropies on 
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an angular scale -(9 is determined by the value of 5p/ p on the corresponding length 
scale \{d) at teq- 

T^ ' 3 p 
On scales of clusters we know that (for ^2 = 1 and h = 1/2) 



■(^) = 77-(M^),W- (4-46) 



(^-A (A(^), te,) ^ ^(te,)-^ ^ 10-\ (4.47) 

\ P / CDM 

using the fact that today on cluster scales Sp/p ~ 1. The bounds on 6T/T on 
small angular scales are 

^(^)10-\ (4.48) 

consistent with the predictions for a CDM model, but inconsistent with those of a 
Q = Qb = 1 model, according to which we would expect anisotropics of the order 
of 10^^. This is yet another argument in support of the existence of nonbaryonic 
dark matter. 

To conclude this subsection, let us briefly discuss Newtonian perturbations 
during the radiation-dominated epoch. We consider matter fluctuations with Cg = 
in a smooth relativistic background. In this case. Equation (4.32) becomes 

5k + 2H6k-A7TGpn,6k = 0, (4.49) 

where pm denotes the average matter energy density. The Hubble expansion pa- 
rameter obeys 

H' = ^{Pm+Pr), (4.50) 

with pr the background radiation energy density. For t < teq, pm is negligible in 
both (4.49) and (4.50), and (4.49) has the general solution 

4(t) = cilogt + C2. (4.51) 

In particular, this result implies that CDM perturbations which enter the Hubble 
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radius before teq have an amplitude which grows only logarithmically in time until 
teq- This is sometimes called the Meszaros effect. 

4.4. Relativistic Theory of Cosmological Perturbations 

On scales larger than the Hubble radius (A > t) the Newtonian theory of cos- 
mological perturbations obviously is inapplicable, and a general relativistic analysis 
is needed. On these scales, matter is essentially frozen in comoving coordinates. 
However, space-time fluctuations can still increase in amplitude. 

In principle, it is straightforward to work out the general relativistic theory of 
linear fluctuations^^-*. We linearize the Einstein equations 

G^^ = 8nGT^^ (4.52) 

(where G^i, is the Einstein tensor associated with the space-time metric g^^, and 
T^jy is the energy-momentum tensor of matter) about an expanding FRW back- 
ground {glJ, (f^^^): 

9^^u{x, t) = gfjit) + Vfe ^) 

^(x,t) = i^^^\t) + 5i^{x,t) 

and pick out the terms linear in h^i, and 5ip to obtain 

5G^^ = SttGST^,. (4.54) 

In the above, h^i, is the perturbation in the metric and Sip is the fluctuation of the 
matter fleld (f. We have denoted all matter flelds collectively by (/?. 

In practice, there are many complications which make this analysis highly 
nontrivial. The flrst problem is "gauge invariance" ^^-^ Imagine starting with a ho- 
mogeneous FRW cosmology and introducing new coordinates which mix x and t. 
In terms of the new coordinates, the metric now looks inhomogeneous. The inho- 
mogeneous piece of the metric, however, must be a pure coordinate (or "gauge") 
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Figure 22: Sketch of how two choices of the mapping from the background space- 
time manifold TVIq to the physical manifold Ai induce two different coordinate 
systems on Ai. 

artefact. Thus, when analyzing relativistic perturbations, care must be taken to 
factor out effects due to coordinate transformations. 

The issue of gauge dependence is illustrated in Fig. 22. A coordinate system on 
the physical inhomogeneous space-time manifold Ai can be viewed as a mapping V 
of an unperturbed space-time A^o iiito Ai. A physical quantity Q is a geometrical 
function defined on TW. There is a corresponding physical quantity ^'^'Q defined 
on A^o- In the coordinate system given by V, the perturbation 6Q of Q at the 
space-time point p e A^ is 

SQ{p) = Q{p)-^^'^Q{D-\p)). (4.55) 

However, in a second coordinate system T> the perturbation is given by 

SQ{p) = Q{p)-(^')Q{V-\p)). (4.56) 
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The difference 

AQ{p) = 5Q{p) - 6Q{p) (4.57) 

is obviously a gauge artefact and carries no physical meaning. 

There are various methods of dealing with gauge artefacts. The simplest and 
most physical approach is to focus on gauge invariant variables, i.e., combinations 
of the metric and matter perturbations which are invariant under linear coordinate 
transformations. 

The gauge invariant theory of cosmological perturbations is in principle straight- 
forward, although technically rather tedious. In the following I will summarize the 
main steps and refer the reader to Ref. 11 for the details and further references (see 
also Ref. 81 for a pedagogical introduction and Refs. 82-87 for other approaches). 

We consider perturbations about a spatially flat Friedmann- Robertson- Walker 
metric 

ds^ = a^{r]){d7f -dx^) (4.58) 

where r] is conformal time (related to cosmic time t by a{ri)dr] = dt). A scalar 
metric perturbation (see Ref. 88 for a precise deflnition) can be written in terms 
of four free functions of space and time: 

Scalar metric perturbations are the only perturbations which couple to energy 
density and pressure. 

The next step is to consider infinitesimal coordinate transformations 

x^'' =x^ + i^ (4.60) 

which preserve the scalar nature of Sg^p and to calculate the induced transforma- 
tions of 0, ■?/', B and E. Then we find invariant combinations to linear order. (Note 
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that there are in general no combinations which are invariant to all orders®^).) 
After some algebra, it follows that 

^ = (f) + a-^[{B-E')a]' 



a 



^ = W, (B-E') 

a 



(AM] 



are two invariant combinations. In the above, a prime denotes differentiation with 
respect to tj. 

There are various methods to derive the equations of motion for gauge invari- 
ant variables. Perhaps the simplest way^^-* is to consider the linearized Einstein 
equations (4.54) and to write them out in the longitudinal gauge defined by 

B = E = (4.62) 

and in which $ = and ^ = ■?/', to directly obtain gauge invariant equations. 

For several types of matter, in particular for scalar field matter, the perturba- 
tion of Tpjy has the special property 

5Ti ~ 5) (4.63) 

which imples $ = \E'. Hence, the scalar-type cosmological perturbations can in 
this case be described by a single gauge invariant variable. The equation of motion 
takes the form^*^'^'^*^-* 



where 

2 //-!$ + $ 

The variable w = p/p (with p and p background pressure and energy density 
respectively) is a measure of the background equation of state. In particular, on 
scales larger than the Hubble radius, the right hand side of (4.64) is negligible, and 
hence ^ is constant. 

53 



The result that .^ = is a very powerful one. Let us first imagine that the 
equation of state of matter is constant, i.e., w = const. In this case, ^ = implies 

$(t) = const , (4.66) 

i.e., this gauge invariant measure of perturbations remains constant outside the 
Hubble radius. 

Next, consider the evolution of $ during a phase transition from an initial 
phase with w = Wi to a phase with w = Wf. Long before and after the transition, 
$ is constant because of (4.66), and hence ^ = becomes 

+ $ = const , (4.67) 



1 + w 



In order to make contact with matter perturbations and Newtonian intuition, it 
is important to remark that, as a consequence of the Einstein constraint equations, 
at Hubble radius crossing $ is a measure of the fractional density fiuctuations: 

^k,tHik))r^^-^{k,tHik)). (4.68) 

(Note that the latter quantity is approximately gauge invariant on scales smaller 
than the Hubble radius). 
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5. Inflationary Universe Scenarios 

5.1. Preliminaries 

Cosmological inflation"^^-' is a period in time during which the Universe is ex- 
panding exponentially, i.e., 

a{t) = e^^ (5.1) 

with constant Hubble expansion rate H. From the FRW equations (2.7) and (2.9) 
it follows that the condition for inflation (in the context of Einstein gravity in a 
spatially fiat Universe) is an equation of state for matter with 

P=-P, (5.2) 

which neccessitates abandoning a description of matter in terms of an ideal gas. 

As was indicated in Section 3.1, it is possible to achieve inflation if matter is 
described in terms of scalar fields, provided that at some period 

V^' « vi^) 

(5.3) 
(Vv^)2 « V{v) 

(see Eqs. (3.9) which give the equation of state for scalar field matter). 



Figure 23: A sketch of two potentials which can give rise to inflation. 
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Two examples which can give inflation are shown in Fig. 23. In (a), inflation 
occurs at the stable flxed point ^p{x,ti) = = 'p{x,ti). However, this model is 
ruled out by observation: the inflationary phase has no ending. V^(0) acts as a 
permanent nonvanishing cosmological constant. In (b), a flnite period of inflation 
can arise if ip{x) is trapped at the local minimum (/? = with ^{x) = 0. However, 
in this case ip{x) can make a sudden transition at some time tji > ti through 
the potential barrier and move to ^p{x) = a. Thus, for ti < t < t^ the Universe 
expands exponentially, whereas for t > tji the contribution of if to the expansion 
of the Universe vanishes and we get the usual FRW cosmology. There are three 
obvious questions: why does the fleld start out at (/? = 0, how does the transition 
occur and why should the scalar fleld have V{ip) = at the global minimum? In 
the following section the flrst two questions will be addressed. The third question is 
part of the cosmological constant problem for which there is as yet no convincing 
explanation. Before studying the dynamics of the phase transition, we need to 
digress and discuss flnite temperature effects. 

5.2. Finite Temperature Field Theory 

The evolution of particles in vacuum and in a thermal bath are very different. 
Similarly, the evolution of flelds changes when coupled to a thermal bath. Under 
certain conditions, the changes may be absorbed in a temperature dependent po- 
tential, the flnite temperature effective potential. Here, a heuristic derivation of 
this potential will be given. The reader is referred to Ref. 8 or to the original 
articles^^-* for the actual derivation. 

We assume that the scalar fleld v^(x, t) is coupled to a thermal bath which is 
represented by a second scalar fleld ^(x, t) in thermal equilibrium. The Lagrangian 
for {f is 

C = Id^^d^v - V{v) - \\fy^ , (5.4) 

where A is a coupling constant. The action from which the equations of motion 
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are derived is 

S = I (fx ^^C (5.5) 



where g is the determinant of the metric (2.2). The resulting equation of motion 
for V5(x, t) is 

(^ + 3/70 -a-2 vV = -'^''M-AV^V- (5.6) 

If il) is in thermal equilibrium, we may replace ilP" by its thermal expectation value 
< ■i/'^ >y_ Now, 

< ^2 >y^ t2 (5.7) 

which can be seen as follows: in thermal equilibrium, the energy density of ip equals 
that of one degree of freedom in the thermal bath. In particular, the potential 
energy density V{ip) of ip is of that order of magnitude. Let 

l^(^) = A^V-' (5.8) 

with a coupling constant A^ which we take to be of the order 1 (if A^ is too 
small, if) will not be in thermal equilibrium). Since the thermal energy density is 
proportional to T^, (5.7) follows. (5.6) can be rewritten as 



i) + "iH^p - a"2 y2 ^ = -\4(v9), (5.9) 

where 

Vt{v) = V{v) + ^ATV (5.10) 

is called the finite temperature effective potential. Note that in (5.10), A has been 
rescaled to absorb the constant of proportionality in (5.7). 
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These considerations will now be applied to Example A, a scalar field model 
with potential 

V{^) = \\{^^-a^f (5.11) 

(o" is called the scale of symmetry breaking). The finite temperature effective 
potential becomes (see Fig. 24) 

For very high temperatures, the effective mass term is positive and hence the 
energetically favorable state is < (/? >= 0. For very low temperatures, on the other 
hand, the mass term has a negative sign which leads to spontaneous symmetry 
breaking. The temperature at which the mass term vanishes defines the critical 
temperature Tg 

Te = A-V2aV2^. (5.13) 



Figure 24: The finite temperature effective potential for Example A. 



As Example B, consider a theory with potential 

with ^a > 6 > 0. The finite temperature effective potential is obtained by adding 
gAT^y?^ to the right hand side of (5.14). Vr(v') is sketched in Fig. 25 for various 
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values of T. The critical temperature Tc is defined as the temperature when the 
two minima of Vt{^p) become degenerate. 



Figure 25: The finite temperature effective potential for Example B. 

It is important to note that the use of finite temperature effective potential 
methods is only legitimate if the system is in thermal equilibrium. This point 
was stressed in Refs. 92 and 93, although the conclusion should be obvious from 
the derivation given above. To be more precise, we require the ip field to be in 
thermal equilibrium and the coupling constant A of (5.4) which mediates the energy 
exchange between the (f and ip fields to be large. However, as shown in Chapter 
4, in inflationary Universe models, the observational constraints stemming from 
the amplitude of the primordial energy density fluctuation spectrum force the self 
coupling constant A of (/? to be extremely small. Since at one loop order, the 
interaction term ■^X^p^ilP' induces contributions to A, it is unnatural to have A 
very small and A unsuppressed. Hence, in many inflationary Universe models - 
in particular in new infiation^'^-' and in chaotic infiation^^-* - finite temperature 
effective potential methods are inapplicable. 
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5.3. Phase Transitions 

The temperature dependence of the finite temperature effective potential in 
quantum field theory leads to phase transitions in the very early Universe. These 
transitions are either first or second order. 

Example A of the previous section provides a model in which the transition is 
second order (see Fig. 24). For T ^ T^, the expectation value of the scalar field ip 
vanishes at all points x in space: 

<V?(x)>=0. (5.15) 

For T < Tc, this value of < v'(x) > becomes unstable and < v'(x) > evolves 
smoothly in time to a new value ±0". The direction is determined by thermal and 
quantum fiuctuations and is therefore not uniform in space. There will be domains 
of average radius C,{t) in which < ip{x) > is coherent. By causality, the coherence 
length is bounded from above by the horizon. However, typical values of C,{t) are 
proportional to X^^a^^ if cp was in thermal equilibrium before the phase transition. 

In condensed matter physics, a transition of the above type is said to proceed 
by spinodal decomposition^^-* , triggered by a rapid quench. 

In Example B of the previous section, (see Fig. 25) the phase transition is first 
order. For T > Tc, the expectation value < v^(a;) > is approximately 0, the min- 
imum of the high temperature effective potential. Provided the zero temperature 
potential has a sufficiently high barrier separating the metastable state ip = from 
the global minimum (compared to the energy density in thermal fiuctuations at 
T = Tc), then v'(x) will remain trapped at (^ = also for T < Tc. In the notation 
of Ref. 96, the field ip is trapped in the false vacuum. After some time (determined 
again by the potential barrier), the false vacuum will decay by quantum tunnelling. 

Tunnelling in quantum field theory was discussed in Refs. 96-99 (for reviews 
see e.g., Refs. 100 and 8). The transition proceeds by bubble nucleation. There is 
a probability per unit time and volume that at a point x in space a bubble of "true 
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vacuum" v^(x) = a will nucleate. The nucleation radius is microscopical. As long 
as the potential barrier is large, the bubble radius will increase with the speed of 
light after nucleation. Thus, a bubble of (^ = a expands in a surrounding "sea" of 
false vacuum (/? = 0. 

To conclude, let us stress the most important differences between the two types 
of phase transitions discussed above. In a second order transition, the dynamics 
is determined mainly by classical physics. The transition occurs homogeneously 
in space (apart from the phase boundaries which - as discussed below - become 
topological defects), and < ^{x) > evolves continuously in time. In first order 
transitions, quantum mechanics is essential. The process is extremely inhomoge- 
neous, and < ^p{x) > is discontinuous as a function of time. As we shall see in 
the following sections, the above two types of transitions are the basis of various 
classes of inflationary Universe models. 

5.4. Models of Inflation 

At this stage we have established the formalism to be able to discuss models of 
inflation. I will focus on "old inflation," "new inflation"" and "chaotic inflation." 
There are many other attempts at producing an inflationary scenario, but there is 
as of now no convincing realization. 

Old Inflation 

The old inflationary Universe model^^'^^^-* is based on a scalar fleld theory which 
undergoes a flrst order phase transition. As a toy model, consider a scalar fleld 
theory with the potential V{ip) of Example B (see Fig. 25). Note that this potential 
is fairly general apart from the requirement that V{a) = 0, where (/? = a is the 
global minimum of V{ip). This condition is required to avoid a large cosmological 
constant today (no inflationary Universe model manages to circumvent or solve 
the cosmological constant problem). 

For fairly general initial conditions, v^(x) is trapped in the metastable state 
(/9 = as the Universe cools below the critical temperature Tc. As the Universe 
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expands further, all contributions to the energy-momentum tensor T^i, except for 
the contribution 

T^u ~ V{(p)g^i, (5.16) 

redshift. Hence, the equation of state approaches p = —p, and inflation sets in. 
Inflation lasts until the false vacuum decays. During inflation, the Hubble constant 
is given by 

H^ = ^V{0). (5.17) 



Figure 26: A sketch of the spatially inhomogeneous distribution of ip in Old 
Inflation. 

After a period T^^, where T is the tunnelling decay rate, bubbles oi cp = a 
begin to nucleate in a sea of false vacuum ip = 0. For a sketch of the resulting 
inhomogeneous distribution of ipix) see Fig. 26. Note that inflation stops after 
bubble nucleation. 

The time evolution in old inflation is summarized in Fig. 27. We denote the 
beginning of inflation by tj (here ti ~ tc), the end by tji (here Ir ^ tc + F^^). 

It was immediately realized that old inflation has a serious "graceful exit" 
problem^^^-* . The bubbles nucleate after inflation with radius r <^ 2t^ and would 
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Figure 27: Phases in tlie old and new inflationary Universe. 

today be much smaller than our apparent horizon. Thus, unless bubbles percolate, 
the model predicts extremely large inhomogeneities inside the Hubble radius, in 
contradiction with the observed isotropy of the microwave background radiation. 

For bubbles to percolate, a sufliciently large number must be produced so that 
they collide and homogenize over a scale larger than the present Hubble radius. 
However, with exponential expansion, the volume between bubbles expands expo- 
nentially whereas the volume inside bubbles expands only with a low power. This 
prevents percolation. 

New Inflation 

Because of the graceful exit problem, old inflation never was considered to be a 
viable cosmological model. However, soon after the seminal paper by Guth, Linde 
and independently Albrecht and Steinhardt put forwards a modifled scenario, the 
New Inflationary Universe^^-* (see also Ref. 103). 

The starting point is a scalar fleld theory with a double well potential which 
undergoes a second order phase transition (Fig. 24). V{ip) is symmetric and ip = 
is a local maximum of the zero temperature potential. Once again, it was argued 
that flnite temperature effects conflne v^(x) to values near ip = aX temperatures 
T > Tc- For T < Tg, thermal fluctuations trigger the instability of v^(x) = and 
(/9(x) evolves towards (/? = ±a by the classical equation of motion 

'P + 3Hip-a-^\y^ip = -V'{ip). (5.18) 
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Figure 28: A sketch of the spatial distribution of ip in New Inflation after the 
transition. The symbols + and — indicate regions where ip = +a and if = —a 
respectively. 

The transition proceeds by spinodal decomposition (see Fig. 28) and hence 
V?(x) will be homogeneous within a correlation length. The analysis will be confined 
to such a small region. Hence, in Eq. (5.18) we can neglect the spatial gradient 
terms. Then, from (3.9) we can read off the induced equation of state. The 
condition for inflation is 



02«\/(v9), (5.19) 



i.e. slow rolling. 



Often, the "slow rolling" approximation is made to find solutions of (5.18). 
This consists of dropping the ip term. In this case, (5.18) becomes 

3Hip =-V'{ip). (5.20) 

As an example, consider a potential which for |(/9| < o" has the following expansion 
near ip = 

Vip) = Vo- ^m'^ip'^ . (5.21) 
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With the above V{(f)), (5.20) has the solution 

(171 \ 
3^V ^^'^^^ 

(taking H = const which is a good approximation). Thus, provided m <^ \^H , (p 
is indeed smaller than the other terms in (5.18) and the slow rolling approximation 
seems to be satisfied. 

However, the above conclusion is premature^^'*-' . Equation (5.18) has a second 
solution. For m > H the solution is 

^(t) ~ v^(0)e"^* (5.23) 

and dominates over the previous one. This example shows that the slow rolling 
approximation must be used with caution. Here, however, the conclusion remains 
that provided m <^ H, then the model produces enough inflation to solve the 
cosmological problems. 

There is no graceful exit problem in the new inflationary Universe. Since the 
spinodal decomposition domains are established before the onset of inflation, any 
boundary walls will be inflated outside the present Hubble radius. 

The condition m^ <^ 3H^ which must be imposed in order to obtain inflation, 
is a flne tuning of the particle physics model - the flrst sign of problems with this 
scenario. Consider e.g. the model (5.11). By expanding V{ip) about (/? = we can 
determine both H and m in terms of A and a. In order that m^ < 3H^ be satisfled 
we need 

/ 1x1/2 

c^> (^1 rripi, (5.24) 

which is certainly an unnatural constraint for models motivated by particle physics. 

Let us, for the moment, return to the general features of the new inflationary 
Universe scenario. At the time tc of the phase transition, ip{t) will start to move 
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from near (/? = towards either ±a as described by the classical equation of 
motion, i.e. (5.22). At or soon after tc, the energy-momentum tensor of the 
Universe will start to be dominated by V{ip), and inflation will commence, ti shall 
denote the time of the onset of inflation. Eventually, 0(t) will reach large values for 
which (5.21) is no longer a good approximation to V{ip) and for which nonlinear 
effects become important. The time at which this occurs is ts- For t > ts , <^{t) 
rapidly accelerates, reaches ±(T, overshoots and starts oscillating about the global 
minimum of V{ip). The amplitude of this oscillation is damped by the expansion of 
the Universe and (predominantly) by the coupling of ip to other fields. At time t^i, 
the energy in (f drops below the energy of the thermal bath of particles produced 
during the period of oscillation. 

The evolution of ip{t) is sketched in Fig. 29. The time period between ts 
and tji is called the reheating period and is usually short compared to the Hubble 
expansion time. The time evolution of the temperature T of the thermal radiation 
bath is also shown in Fig. 29. 



Figure 29: Evolution of ip(t) and T(t) in the new inflationary Universe. 



Reheating in inflationary Universe models has been considered in Refs. 105- 
108. One way to view the process is as follows '•*. Consider a second scalar 
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field ip coupled to (/? via the interaction Lagrangian 

Ci = ^gv^^^ . (5.25) 

Then, an oscillating (p{t) will act as a time dependent mass with periodic variations 
in the equation of motion for ip 

4 + 3H^k + {ml + k'^a^^t) + gy^\t)) ^k = (5.26) 

where we have neglected nonlinear terms and expanded ijj into Fourier modes ipk- If 
the expansion of the Universe can be neglected and for periodic v^^(t), the above is 
the well known Mathieu equation^^^-* whose solutions have instabilities for certain 
values of k. These instabilities correspond to the production oiip particles with well 
determined momenta '-*. These particles eventually equilibrate and regenerate 
a thermal bath. 

For t > t^, the Universe is again radiation dominated. Hence, the stages of 
the new inflationary Universe are the same as for old inflation (Fig. 27). There 
is a useful order of magnitude relation between the scale of symmetry breaking a 
and H. From 

H'=^-^ViO) (5.27) 

and from the form of the potential (see (5.11)) it follows that 

^"l ~ aV2 (^) . (5.28) 

'^ J \mpij 

In particular, for a ~ lO^^GeV (typical scale of grand unification) and A ~ 1 we 
obtain H ~ lOi^GeV. 

The new infiationary Universe model - although it was for a long time presented 
as a viable model - suffers from severe fine tuning and initial condition problems. 
In (5.24) we encountered the first of these problems: in order to obtain enough 
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inflation, the potential must be fairly flat near ip = 0. A more severe problem 
will be derived in Chapter 4: Inflationary Universe models generate energy density 
perturbations. The steeper the potential, the larger the density perturbations. For 
a potential which near ip = ot ip = H has the following expansion 

V{ip) = V{0)-Xip'^, (5.29) 

the density perturbations conflict with observations unless (see later) 



This in itself is an unexplained small number problem. However, even if we 
were willing to accept this we would run into initial condition problems^^'^"^-* . For 
the new inflationary Universe to proceed in the way outlined above, it is essential 
that the field (f be in thermal equilibrium with other fields. This implies that the 
constant g coupling (/? to other fields should not be too small. However, a coupling 
term of the form (5.25) induces one loop quantum corrections to the self coupling 
constant A of the order g'^. Hence, the constraint A < 10^^^ implies a constraint 
g < 10^^. Thus, (p will not be in thermal equilibrium at tc, and hence there will 
be no thermal forces which localize ip close to (/? = 0. 

Note that the above problem is not an artifact of using quartic potentials 
such as (5.29). Similar constraints would arise in other {e.g. quadratic) models. 
However, (5.29) was long considered to be the prototypical shape of V{(p) for small 
values of cp since it is the shape which arises in Coleman- Weinberg^ ^^•' models. 

In the absence of thermal forces which constrain ip to start close to (/? = 0, the 
only constraints on (p— at least using classical physics alone - come from energetic 
considerations. Obviously, it is unnatural to assume that at the initial time ti the 
energy density in ip exceeds the energy density of one degree of freedom of the 
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thermal bath at time U (temperature Tj). This imphes 

\y ^{ilM' <YqTS\U) (5.31) 

In particular, for the double well potential of (5.11), (5.31) implies that ^^(x, ti) 
can be of the order 

^(x,t,)~A-i/4Ti (5.32) 

which for T^ > o" is much larger than a. In a weakly coupled model, the only 
natural time to impose initial conditions on v'(x) is the Planck time, i.e. Ti ~ irtpi. 
Hence, the initial conditions allow and in fact suggest 



(x,ti)~A ^^^rupi :$> rupi forTi~Tpi, (5.33) 



V 



These observations lead to the chaotic inflation scenario, the only of the original 
inflationary Universe models which can still be considered as a viable scenario 
today. 

Chaotic Inflation 

Chaotic inflation^^-' is based on the observation that for weakly coupled scalar 
fields, initial conditions which follow from classical considerations alone lead to 
very large values of ip{x) (see (5.32)). 

Consider a region in space where at the initial time v^(x) is very large, homo- 
geneous (we will make these assumptions quantitative below) and static. In this 
case, the energy-momentum tensor will be immediately dominated by the large 
potential energy term and induce an equation of state p ~ —p which leads to infla- 
tion. Due to the large Hubble damping term in the scalar field equation of motion, 
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V9(x) will only roll very slowly towards (/? = 0. The kinetic energy contribution 
to T^i, will remain small, the spatial gradient contribution will be exponentially 
suppressed due to the expansion of the Universe, and thus inflation persists. This 
is a brief survey of the chaotic inflation scenario. Note that in contrast to old and 
new inflation, no initial thermal bath is required. Note also that the precise form 
of V{ip) is irrelevant to the mechanism. In particular, V{ip) need not be a double 
well potential. This is a signiflcant advantage, since for scalar flelds other than 
Higgs flelds used for spontaneous symmetry breaking, there is no particle physics 
motivation for assuming a double well potential, and since the inflaton (the fleld 
which gives rise to inflation) cannot be a conventional Higgs fleld due to the severe 
flne tuning constraints. 

Let us consider the chaotic inflation scenario in more detail. For simplicity, 
take the potential 

\/(v9)=^mV (5.34) 

and consider a region in space in which ip{x,ti) is sufficiently homogeneous. To be 
specific, we require 

ia-2(t^)| V^(x,t,)|2 « V{^{x,U)) (5.35) 

over a region of size di 

di>3H-^{ti). (5.36) 

We also require that the kinetic energy be negligible at the initial time tj, 

0(x,t.)2 < V(v^(x,ti)) , (5.37) 

although this assumption can be relaxed without changing the results^^^-*. From 
(5.35) and (5.37) it follows that at ti the equation of state is infiationary, i.e. p{ti) ~ 
—p{ti). Condition (5.36) ensures that no large inhomogeneities can propagate from 
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outside to the center of the region under consideration. With these approximations, 
the equation of motion for if becomes 



(p + "iH^ = -m^ip (5.38) 



with 

H=(^-^X^' —V. (5.39) 



3 J rripi 

Since we expect ^^(x, t) to be changing slowly, we make the slow rolling approx- 
imation 

3/70 = -mV (5.40) 

which gives 

= -(Y2^1 mmpi (5.41) 

and shows that the approximation is self consistent. In order to get inflation, we 
require 

l^^<Un\^ (5.42) 

which (by (5.41)) is satisfied if 

/ 1x1/2 

V^ > ( Y^ J ^pi ■ (5.43) 

In order to obtain a period r > 50iJ^^(tj) of inflation, a slightly stronger condition 
is needed: 

LP > Srupi . (5.44) 

With chaotic inflation, the initial hope that grand unified theories could pro- 
vide the answer to the homogeneity and flatness problems has been abandoned. 
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The inflaton is introduced as a new scalar field (with no particular particle physics 
role) which is very weakly coupled to itself and to other fields (see Ref. 112 for an 
attempt to couple the inflaton to non-grand unifled particle physics). In supergrav- 
ity and in superstring inspired models there are scalar flelds which are candidates 
to be the inflaton. I refer the reader to Refs. 113-115 for a discussion of this issue. 
However, even in such models the time when the inflaton >^ decouples from the 
rest of physics is the Planck time ti = tpi- Thus, the chaotic inflation scenario is 
often called primordial inflation. 

Chaotic inflation is a much more radical departure from standard cosmology 
than old and new inflation. In the latter, the inflationary phase can be viewed 
as a short phase of exponential expansion bounded at both ends by phases of 
radiation domination. In chaotic inflation, a piece of the Universe emerges with an 
inflationary equation of state immediately after the quantum gravity epoch. 

The chaotic inflationary Universe scenario has been developed in great detail 
(see e.g., Ref. 116 for a recent review). One important addition is the inclusion of 
stochastic noise^^^-* in the equation of motion for ip in order to take into account 
the effects of quantum fluctuations. It can in fact be shown that for sufficiently 
large values of 1^9 1, the stochastic force terms are more important than the classical 
relaxation force V'{ip). There is equal probability for the quantum fluctuations to 
lead to an increase or decrease of \ip\. Hence, in a substantial fraction of comoving 
volume, the fleld ip will climb up the potential. This leads to the conclusion that 
chaotic inflation is eternal. At all times, a large fraction of the physical space 
will be inflating. Another consequence of including stochastic terms is that on 
large scales (much larger than the present Hubble radius), the Universe will look 
extremely inhomogeneous. 

General Comments 

Old, new and choatic inflation are all based on the use of new fundamental 
scalar flelds which cannot be the Higgs fleld of an ordinary gauge theory. Instead 
of introducing new physics via scalar flelds - an approach which makes the cos- 
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mological constant problem more severe - it is possible to look for realizations of 
inflation based on some alternative new physics which do not invoke fundamental 
scalar fields. 

One possibility is to consider modifications of Einstein gravity which can lead 
to inflation. In fact, the first model of infiation^^-* was based on considering an 
action for gravity of the form 



S= (rx^/^{R + eR^), (5.45) 



where R is the Ricci scalar of the metric g^^. As shown first by Whitt^^®-*, the 
equations of motion resulting from this action are the same as those from Einstein 
gravity in the presence of a scalar field with a special potential which allows for 
chaotic inflation. The relationship is obtained via a conformal transformation. 

Since perturbative quantum gravity calculations and considerations based on 
quantum field theory in curved space-time all point to the presence of higher deriva- 
tive terms in the action for gravity, it is not unlikely that successful inflation will 
come from the gravity sector. A recently proposed theory -*, in which the curva- 
ture is bounded for all solutions, predicts that our Universe will have started out 
in an inflationary period. 

Another interesting possibility is that inflation is a result of new physics as- 
sociated with a unified theory of all forces such as superstring theory. Interesting 
speculations along these lines have recently been made in Ref. 119. 

As has hopefully become clear, inflation is a nice idea which solves many prob- 
lems of standard big bang cosmology. However, no convincing realization of infla- 
tion which does not involve unexplained small numbers has emerged (for a general 
discussion of this point see Ref. 120). 

It is important to distinguish between models of inflation which are self con- 
sistent and those which are not. We have shown that new inflation is not self 
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consistent, whereas chaotic inflation is. One of the key issues involves initial con- 
ditions. In new inflation, the initial conditions required can only be obtained if 
the inflaton field is in thermal equilibrium above the critical temperature, which 
however is not possible because of the density fluctuation constraints on coupling 
constants. 

In chaotic inflation, it can be shown that - provided the spatial sections are 
flat - a large phase space of initial conditions (much larger than is apparent from 
(5.35) and (5.37)) gives chaotic inflation^^^^^^^-*, whereas the probability to relax 
dynamically^^^-* to fleld conflgurations which give new inflation (this possibility is 
only available in double well potentials) is negligibly small. 

5.5. Generation and Evolution of Fluctuations 

Preliminaries 

In this section, the origin of the primordial density perturbations required to 
seed galaxies will be discussed within the context of inflationary Universe models. 
From Chapter 3, we recall the basic reason why in inflationary cosmology a causal 
generation mechanism is possible: comoving scales of cosmological interest today 
originate inside the Hubble radius early in the de Sitter period (see Fig. 11). 
Hence, it is in principle possible that density perturbations on these scales can be 
generated by a causal mechanism at very early times. 

First, let us demonstrate why the Hubble radius H^^(t) is the length scale of 
relevance in these considerations. Consider a scalar fleld theory with action 



\d^^d^^ - v{^) 



(5.46) 



S{ip) = d x^f^ 

The resulting equation of motion is 

(p + 'iE(p - a-^vV = -V\^) . (5.47) 

The second term on the left hand side is the Hubble damping term, the third repre- 
sents microphysics (spatial gradients). To simplify the consideration, assume that 
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V'{lp) = 0. Then, the time evolution is influenced by microphysics and gravity. For 
plane wave perturbations with wave number k, the gravitational force is propor- 
tional to H^cp whereas the microphysical force is a~'^k^Lp. Thus, for ak^^ < H^^, 
i.e., on length scales smaller than the Hubble radius, microphysics dominates, 
whereas for ak^^ > H~^, i.e., on length scales larger than the Hubble radius, the 
gravitational drag dominates. 

Based on the above analysis we can formulate the main idea of the fluctuation 
analysis in inflationary cosmology. In linear order, all Fourier modes decouple. 
Hence, we fix a mode with wave number k. There are two very different time 
intervals to consider. Let ti{k) be the time when the scale crosses the Hubble radius 
in the de Sitter phase, and tf{k) the time when it reenters the Hubble radius after 
inflation (see Fig. 30). The first period lasts until t = ti{k). In this time interval 
microphysics dominates. We shall demonstrate that quantum fluctuations generate 
perturbations during this period"^^""^®'^^^-*. The second time interval is ti{k) < t < 
tf{k). Now microphysics is unimportant and the evolution of perturbations is 
determined by gravity. At ti{k), decoherence sets in; the quantum mechanical 
wave functional evolves as a statistical ensemble of classical configurations after 
this time. 

It is possible to give a heuristic derivation of the shape of the spectrum of 
cosmological perturbations in an inflationary Universe based on simple geometri- 
cal arguments^^-* . Consider flrst the process of generation of fluctuations. If the 
generation mechanism produces inhomogeneities at all times on a flxed physical 
wavelength with an amplitude determined by the "Hawking temperature" Tjj of 
de Sitter space^^^-*, 

Th = ^, (5.48) 

then the evolution of perturbations on different scales between when they are 
formed and when they leave the Hubble radius at times ti{k) is related by time 
translation, and the amplitude of the fluctuations when measured at time ti{k) will 
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Figure 30: Sketch (physical distance Xp versus time t) of the evolution of two 
fixed comoving scales labelled ki and k2 in the inflationary Universe. 



be independent of k: 

SM 



{k,ti{k)) =Ai= const. (5.49) 



M 

In fact, the amplitude Ai will be given by the thermal energy associated with 
the Hawking temperature (5.49) divided by the "false vacuum" energy density po 
responsible for inflation 

T-l4 

A ~ — • (5.50) 

PO 

In Chapter 4 it was shown (see (4.68)) that the gauge invariant measure of 
density fluctuations evolves trivially on scales outside of the Hubble radius: it 
changes by a factor which only depends on the equations of state in the inflationary 
and past-inflationary phases. In particular, this factor is independent of k. Hence, 
using (4.68), we arrive at the conclusion that perturbations are independent of 
scale when they enter the Hubble radius 

-—{k,tf{k)) =Af= const. (5.51) 

Hence, we have demonstrated that a scale invariant spectrum of density perturba- 
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tions is a generic feature of an inflationary Universe scenario. 
From (4.67) and (4.68) if follows that 



I+Wjtfjk)) 



On scales which enter the Hubble radius after tgq, 

l+w{tf{k)) = l. (5.53) 

The initial value oi w = p/p during inflation can be determined by using the 
expressions (3.9) for energy density and pressure of the scalar field responsible for 
inflation. The result is 



!+„(,(*))= ».^. (5^54) 

p{U{k)) po 



where the last proportionality follows by dimensional analysis. Thus, combining 

(5.50-5.54) we obtain 

6M 



M 



{k,tf{k))r^l, (5.55) 



which is at least four orders of magnitude too large to conform with the constraints 
from CMB anisotropy measurements. 

The above problem is known as the "fluctuation problem" 127,90) ^^^^ jg common 
to most inflationary Universe models. The only known solutions involve small 
numbers introduced into the particle physics sector. This defeats one of the aims 
of inflation which is to avoid the need for unnaturally small constants. To study this 
problem in more detail we must turn to a quantitative analysis of the generation 
and evolution of fluctuations. 

77 



Quantum Generation of Fluctuations 

The question of the origin of classical density perturbations from quantum 
fluctuations in the de Sitter phase of an inflationary Universe is a rather subtle 
issue. Starting from a homogeneous quantum state {e.g., the vacuum state in the 
FRW coordinate frame at time ti, the beginning of inflation), a naive semiclassical 
anaylsis would predict the absence of fluctuations since 

< i/j\Tf^u{x)\i{j >= const . (5.56) 

However, as a simple thought experiment shows, such a naive analysis is in- 
appropriate. Imagine a local gravitational mass detector D positioned close to a 
large mass M which is suspended from a pole (see Fig. 31). The decay of an alpha 
particle will sever the cord (at point T) by which the mass is held to the pole and 
the mass will drop. According to the semiclassical prescription 

G^, = 8nG<^\T^,\^P>, (5.57) 

the metric {i.e., the mass measured) will slowly decrease. This is obviously not 
what happens. The mass detector shows a signal which corresponds to one of the 
classical trajectories which make up the state {ip >, a trajectory corresponding to 
a sudden drop in the gravitational force measured. 



Figure 31: Sketch of the thought experiment discussed in the text. 
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The origin of classical density perturbations as a consequence of quantum fluc- 
tuations in a homogeneous state {ip > can be analyzed along similar lines. The 
quantum to classical transition is picking out^^®"^"^^-' one of the typical classical 
trajectories which make up the wave function of \4' >■ We can implement^^^'^'^^-' 
the procedure as follows: Define a classical scalar field 

Vcl{x, t) = ipo{t) + 5ip{x, t) (5.58) 



with vanishing spatial average of 5ip. The induced classical energy momentum 

pel. 



tensor Trj^(x,t) which is the source for the metric is given by 



T;t(ai,t)=V(^rffet)), (5.59) 

where T^i, {'^cl{^,t)) is defined as the canonical energy- momentum tensor of the 
classical scalar field '^dilL^'t)- Unless 5ip vanishes, Tj^, is inhomogeneous. 

For applications to chaotic inflation, we take {ip > to he a Gaussian state with 
mean value v^o(^) 

<^l<^'fet)|^>=V^[](t). (5.60) 

Its width is taken to be the width of the vacuum state of the free scalar field theory 
with mass determined by the curvature of V^(v^) at ipQ. This state is used to define 
the Fourier transform 6ip{k,t) by 

\S^{k)f=<^\\^{k)\^\^> . (5.61) 

The amplitude of 6ip{k) is identified with the width of the ground state wave 
function of the harmonic osciallator (p{k). (Recall that each Fourier mode of a 
free scalar field is a harmonic oscillator). Note that no divergences arise in the 
above construction. In principle, quantum fluctuations contribute a term to ^po{t), 
but this term will be negligible compared to the expectation value (5.60), at least 

79 



in chaotic inflation. The quantum corrections to (5.60) are divergent and must 
be regularized and renormahzed (see e.g., Ref. 133). They are the source of the 
stochastic driving forces in stochastic chaotic inflation. 

By hnearizing (5.59) about v^o(^) we obtain the perturbation of the energy- 
momentum tensor. In particular, the energy density fluctuation Sp{k) is given 
by 

6p{k) = ^Q6(p{k) + V\^^^6^{k) . (5.62) 

To obtain the initial amplitude Ai of (5.49), the above is to to evaluated at the 
time ti{k). 

The computation of the spectrum of density perturbations produced in the de 
Sitter phase has been reduced to the evaluation of the expectation value (5.61). 
First, we must specify the state \ip >. (Recall that in non-Minkowski space-times 
there is no uniquely deflned vacuum state of a quantum fleld theory). We pick the 
FRW frame of the pre-inflationary period. In this frame, the number density of 
particles decreases exponentially. Hence we choose {ip > to he the ground state 
in this frame (see Ref. 134 for a discussion of other choices). ■?/'[(^(k), t], the wave 
functional of {ip >, can be calculated explicitly. It is basically the superposition of 
the ground state wave functions for all oscillators 

ij[^{k),t] = iVexp |-i(27r)"V(t) f d^kuj{k,t)\^{k)\A . (5.63) 

A^ is a normalization constant and uj{k,t) ~ if at t = ti{k). Hence 

(5^(k, t) = {2nf/^a'^/^uj{k, t)-i/2 ~ {2nf/\-^/'^H , t = ti{k) . (5.64) 

Evolution of Fluctuations 

Given the above determination of the intitial amplitude of density perturba- 
tions at the time when they leave the Hubble radius during the de Sitter phase, 
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and the general relativistic analysis of the evolution of fluctuations discussed in 
Section 4.4, it is easy to evaluate the r.m.s. inhomogeneities when they reenter the 
Hubble radius at time tf{k). 

First, we combine (5.62), (5.64), (5.49), (4.10) and (4.15) to obtain 



M J V Po / V Po 

and thus 

A,~nMiM£. (5.66) 

Po 

If the background scalar field is rolling slowly, then 

V'{ipo{ti{m = 3H\ipo{ti{k))\. (5.67) 

Combining (5.66) and (5.67) with (5.52) and (5.54) we get 

SM,^ ,,,, , 3H'^\iPo(ti(k))\ 3H^ 3H^ ,, ^„, 

This result can now be evaluated for specific models of inflation to find the condi- 
tions on the particle physics parameters which give a value 

Af ~ 10~^ (5.69) 

which is required if quantum fluctuations from inflation are to provide the seeds 
for galaxy formation and agree with the CMB anisotropy limits. 

For chaotic inflation with a potential 

\/(v9) = imV, (5.70) 

the dynamics of (/? was analyzed in Section 5b (see in particular (5.43) and (5.44)). 
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We have 

Lp{ti{k)) r^ rupi (5.71) 

and hence 

H{tiik))r~.m-^mip{tiik))r^m. (5.72) 

Therefore, 

^-§ik, tf{k)) ~ 3 f ~ 10— (5.73) 

which imphes that 

m-lO^^QgV (5.74) 

to agree with (5.70). 

Similarly, for a potential of the form 



V{v) = ^A/ (5.75) 



we obtain 



which requires 



^(A:,t;(A;))~10-AV2 (5.76) 



A < 10"^2_ (5 77) 



in order not to conflict with observations. 

The conditions (5.74) and (5.77) require the presence of small parameters in the 
particle physics model. It has been shown quite generally^^^-* that small parameters 
are required if inflation is to solve the fluctuation problem. 
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Discussion 

Let us first summarize the main results of tfie analysis of density fluctuations 
in inflationary cosmology: 

— Quantum vacuum fluctuations in the de Sitter phase of an inflationary Uni- 
verse are the source of perturbations. 

— The quantum perturbations decohere on scales outside the Hubble radius 
and can hence be treated classically. 

— The classical evolution outside the Hubble radius produces a large amplifica- 
tion of the perturbations. In fact, unless the particle physics model contains 
very small coupling constants, the predicted fluctuations are in excess of 
those allowed by the bounds on cosmic microwave anisotropies. 

— Inflationary Universe models generically produce a scale invariant 
Harrison-Zel'dovich spectrum 

-^{k,tf{k)) = const. (5.78) 

It is not hard to construct models which give a different spectrum. All that is 
required is a significant change in H during the period of inflation. 

I have chosen to present the analysis of fluctuations in inflationary cosmology in 
two separate steps in order to highlight the crucial physics issues. Having done this, 
it is possible to step back and construct a unified analysis in which expectation 
values of gauge invariant variables are propagated from t <^ ti{k) to tf{k) in a 
consistent way^"^^'^^-*, and in which the final values of the expectation values of 
quadratic operators are used to construct T^[,(x, t). 

Once inside the Hubble radius, the evolution of the mass perturbations is 
influenced by the damping effects discussed in Section 4.3, which is turn depend 
on the composition of the dark matter. The dominant effects are the Meszaros 
effect and free streaming. 
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On scales which enter the Hubble radius before teg, the perturbations can 
only grow logarithmically in time between tf{k) and tgg- This implies that (up to 
logarithmic corrections), the mass perturbation spectrum is flat for wavelengths 
smaller than X^q, the comoving Hubble radius at teq'- 

— (A,t) ~ const, t<teq,X<Xeq, (5.79) 

whereas on larger scales 

— (A,t)ocA-2. (5.80) 

Equations (5.78) and (5.79) give the power spectrum in an fi = 1 inflationary CDM 
model. 

If the dark matter is hot, then neutrino free streaming cuts off the power 
spectrum at Aj?"^ (see (3.27)). The inflationary HDM and CDM perturbation 
spectra are compared in Fig. 32. 

6. Topological Defects and Structure Formation 

6.1. Classification 

In the previous section we have seen that symmetry breaking phase transitions 
in unifled fleld theories arising in particle physics {e.g., Grand Unifled Theories 
(GUT)^^-*) do not lead, in general, to inflation. In most models, the coupling 
constants which arise in the effective potential for the scalar fleld ip driving the 
phase transition are too large to generate a period of slow rolling which lasts more 
than one Hubble time H^^{t). Nevertheless, there are interesting remnants of the 
phase transition: topological defects. 

Consider a single component real scalar fleld with a typical symmetry breaking 
potential 

V{^) = \x{^^-rff (6.1) 

Unless A ^ 1 there will be no inflation. The phase transition will take place on a 
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Figure 32: The linear tlieory power spectra for inflationary CDM (upper curve) 
and HDM models. The horizontal axis is mass, expressed in units of solar masses. 
Meq is the mass inside the comoving horizon at teq-, and Mps is the mass inside 
the maximal comoving neutrino free streaming volume. 

short time scale r < H^^, and will lead to correlation regions of radius ^ < t inside 
of which if is approximately constant, but outside of which (/? ranges randomly 
over the vacuum manifold M.^ the set of values of (/? which minimizes V{ip) - in our 
example (/? = ir^. The correlation regions are separated by domain walls, regions 
in space where ip leaves the vacuum manifold M. and where, therefore, potential 
energy is localized. Via the usual gravitational force, this energy density can act 
as a seed for structure. 

Topological defects are familiar from solid state and condensed matter sys- 
tems. Crystal defects, for example, form when water freezes or when a metal 
crystallizes •*. Point defects, line defects and planar defects are possible. Defects 
are also common in liquid crystals. They arise in a temperature quench from the 
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disordered to the ordered phase^^^-* . Vortices in "^He are analogs of global cosmic 
strings. Vortices and other defects are also produced during a quench below the 
critical temperature in ^He^"^®-* . Finally, vortex lines also play an important role in 
the theory of superconductivity^^^-* . 

The analogies between defects in particle physics and condensed matter physics 
are quite deep. Defects form for the same reason: the vacuum manifold is topolog- 
ically nontrivial. The arguments which say that in a theory which admits defects, 
such defects will inevitably form, are applicable both in cosmology and in con- 
densed matter physics. Different, however, is the defect dynamics. The motion of 
defects in condensed matter systems is friction-dominated, whereas the defects in 
cosmology obey relativistic equations, second order in time derivatives, since they 
come from a relativistic field theory. 

After these general comments we turn to a classification of topological defects. 
We consider theories with an n-component order parameter ip and with a potential 
energy function (free energy density) of the form (6.1) with 

n 

v'=y: vi ■ (6.2) 

i=l 

There are various types of local and global topological defects (regions of 
trapped energy density) depending on the number n of components of ip. The 
more rigorous mathematical definition refers to the homotopy of Ai. The words 
"local" and "global" refer to whether the symmetry which is broken is a gauge or 
global symmetry. In the case of local symmetries, the topological defects have a 
well defined core outside of which p contains no energy density in spite of nonva- 
nishing gradients Vp: the gauge fields A^ can absorb the gradient, i.e., D^p = 
when dfj,p ^ 0, where the covariant derivative D^ is defined by 

D^ = d^ + ieA^, (6.3) 

e being the gauge coupling constant. Global topological defects, however, have 
long range density fields and forces. 
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Table 1 contains a list of topological defects with their topological character- 
istic. A "v" markes acceptable theories, a "x" theories which are in conflict with 
observations (for r] ~ 10 GeV). 



Theories with domain walls are ruled out^^*^-* since a single domain wall stretch- 
ing across the Universe today would overdose the Universe. Local monopoles are 
also ruled out^^^-* since they would overdose the Universe. Local textures are 
ineffective at producing structures (see Section 6). 

We now describe examples of domain walls, cosmic strings, monopoles and 
textures, focussing on configurations with maximal symmetry. 
Domain walls 

arise in theories with a single real order parameter and free energy density 
given by (6.1). The vacuum manifold of this model consists of two points 

M = {±r]} (6.4) 

and hence has nontrivial zeroth homotopy group: 

no(A^) ^ 1 (6.5) 

(readers not familiar with homotopy groups can simply skip all of the following 
statements involving Iln{Ai). They are not required for an understanding of the 
physics) . 
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To construct a domain wall configuration with planar symmetry (without loss 
of generality the y — z plane can be taken to be the plane of symmetry), assume 
that 

(6.6) 
Lp{x) ~ —7] X <^ —rj" 

By continuity of y?, there must be an intermediate value of x with ^{x) = 0. We 
can take this point to be a; = 0, i.e., 

V9(0)=0. (6.7) 

The set of points with (p = constitute the center of the domain wall. Physically, 
the wall is a thin sheet of trapped energy density. The width w of the sheet is 
given by the balance of potential energy and tension energy. Assuming that the 
spatial gradients are spread out over the thickness w we obtain 

wV{0)=wXr]'^ ^ -7]^ (6.8) 

w 

and thus 

w r^ X-^/^T]-^ . (6.9) 

See Fig. 33 for a sketch of this domain wall configuration. 



Figure 33: A low temperature field configuration containing a domain wall. The 
sketch shows a plane P in space. At positions with a + or — , the value of the 
scalar field is rj or —77, respectively. The solid line is the intersection with P of the 
plane of points in space with (/? = 0. 



A theory with a complex order parameter (n = 2) admits cosmic strings. In 
this case the vacuum manifold of the model is 

M = S\ (6.10) 

which has nonvanishing first homotopy group: 

Ui{M) = Z^l. (6.11) 

A cosmic string is a line of trapped energy density which arises whenever the 
field v^(x) circles Ai along a closed path in space {e.g., along a circle). In this 
case, continuity of (/? implies that there must be a point with (/? = on any sheet 
bounded by the closed path. The points on different sheets connect up to form a 
line overdensity of field energy (see Fig. 34). 



Figure 34: Sketch of the topological argument for the existence of cosmic string 
configurations. Given a field configuration with nontrivial winding along a circle 
normal to the plane of this figure, there must be a point with (/? = on every disk 
bounded by the circle. Three disks are depicted: D, D' and D" , and the respective 
points with (/? = are z, z' and z" . The union of all such points z forms the center 
z{\) of the string. 

To construct a field configuration with a string along the z axis ^ take v^(x) 
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to cover M. along a circle with radius r about the point (x, y) = (0, 0): 

(y9(r,^) ~ r^e*'^, r > r]"^ . (6.12) 

This configuration has winding number 1, i.e., it covers Ai exactly once. Main- 
taining cylindrical symmetry, we can extend (6.12) to arbitrary r 

^(r, ^) = /(r)e^^ (6.13) 

where /(O) = and /(r) tends to rj for large r. The width w can again be found by 
balancing potential and tension energy. The result is identical to the result (6.9) 
for domain walls. 

For local cosmic strings, i.e., strings arising due to the spontaneous breaking 
of a gauge symmetry, the energy density decays exponentially for r ^ r]~^. In this 
case, the energy fi per unit length of a string is finite and depends only on the 
symmetry breaking scale i] 

^i^rf (6.14) 

(independent of the coupling A). The value of /i is the only free parameter in a 
cosmic string model. 

To see how the finiteness of the mass per unit length /i comes about, con- 
sider the simplest theory admitting local strings, the Abelian Higgs model with 
Lagrangean 

£ = ^ D^^D^'^ - Vi^) + \f^,F^- , (6.15) 

where (y? is a complex order parameter with potential (6.1), D^ is the gauge covari- 
ant derivative 



D^ = d^ + ieA^, (6.16) 

the field A^ is a U(l) gauge potential with associated field strength 

F^^ = d^A,-d,,A^, (6.17) 

and e is the gauge coupling constant. 
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For an order parameter configuration (6.12), the gauge fields A^ will take on 
values such that 

D^^ ~ r > 7]-^ (6.18) 

even though d^Lp ^ 0. Hence, the energy density decays exponentially for r ^ ri~^. 
For strings in a global theory (no gauge potential), the spatial gradient energy 
(d^if)'^ cannot be cancelled at large r, and hence the mass per unit length is 
logarithmically divergent as a function of a large r cutoff. 

If the order parameter of the model has three components (n = 3), then 
monopoles result as topological defects. The vacuum manifold is 

M = S^ (6.19) 



and has topology given by 



U2{M)^1. (6.20) 



Given a sphere S is space, it is possible that (/? takes on values in Ai everywhere on 
S, and that it covers Ai once. By continuity, there must be a point in space in the 
interior of S with (/? = 0. This is the center of a point-like defect, the monopole. 

To construct a spherically symmetric monopole with the origin as its center, 
consider a field configuration ip which defines a map from physical space to field 
space such that all spheres Sr in space of radius r ^ i]^^ about the origin are 
mapped onto A4 (see Fig. 35): 

(f : Sr > M 

(6.21) 

(r,^,^) — > (^,V^). 

This configuration defines a winding number one magnitude monopole. 

Domain walls, cosmic strings and monopoles are examples of topological defects . 
A topological defect has a well-defined core, a region in space with ip ^ Ai and 
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hence V{lp) > 0. There is an associated winding number which is quantized, i.e., it 
can take on only integer values. Since the winding number can only change continu- 
ously, it must be conserved, and hence topological defects are stable. Furthermore, 
topological defects exist for theories with global and local symmetries. 

Now, let us consider a theory with a four-component order parameter [i.e., 
n = 4), and a potential given by (6.1). In this case, the vacuum manifold is 

M = S^ (6.22) 

and the associated topology is given by 

Us{M) y^l. (6.23) 

The corresponding defects are called " textures ".^'^'^^-* 



Figure 35: Construction of a monopole: left is physical space, right the vacuum 
manifold. The field configuration maps spheres in space onto Ai. However, a 
core region of space near the origin is mapped onto field values not in Ai. 
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Textures, however, are quite different than the previous topological defects. 
The texture construction will render this manifest (Fig. 36). To construct a 
radially symmetric texture, we give a field configuration ip{x) which maps physical 
space onto Ai. The origin in space (an arbitrary point which will be the center 
of the texture) is mapped onto the north pole N of Ai. Spheres surrounding 
are mapped onto spheres surrounding A^. In particular, some sphere with radius 
rc{t) is mapped onto the equator sphere of Ai. The distance rc{t) can be defined 
as the radius of the texture. Inside this sphere, f{x) covers half the vacuum 
manifold. Finally, the sphere at infinity is mapped onto the south pole oi Ai. The 
configuration v^(x) can be parameterized by^^-* 

V{x,y,z) = cosx(r), sinx(r)-, sinx(r)-, sinx(r)- (6.24) 

in terms of a function x(r) with x(0) = and x(oo) = vr. Note that at all points 
in space, v^(x) lies in At. There is no defect core. All the energy is spatial gradient 
(and possibly kinetic) terms. 

In a cosmological context, there is infinite energy available in an infinite space. 
Hence, it is not necessary that x(^) ^ vr as r -^ oo. We can have 

Xir) -^ Xmax < TT as r ^ OO . (6.25) 

In this case, only a fraction 

nw = t; (6-26) 

of the vacuum manifold is covered: the winding number nw is not quantized. This 
is a reflection of the fact that whereas topologically nontrivial maps from S^ to S^ 
exist, all maps from R^ to S^ can be deformed to the trivial map. 

Textures in R are unstable. For the conflguration described above, the insta- 
bility means that rc{t) — » as t increases: the texture collapses. When rc(t) is 
microscopical, there will be sufficient energy inside the core to cause v^(0) to leave 
At, pass through and equilibrate at x(0) = vr: the texture unwinds. 
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Figure 36: Construction of a global texture: left is physical space, right the 
vacuum manifold. The field configuration is a map from space to the vacuum 
manifold (see text). 

A further difference compared to topological defects: textures are relevant only 
for theories with global symmetry. Since all the energy is in spatial gradients, for a 
local theory the gauge fields can reorient themselves such as to cancel the energy: 

D^^ = 0. (6.27) 

Therefore, it is reasonable to regard textures as an example of a new class of 
defects, semitopological defects. In contrast to topological defects, there is no core, 
and V'(x) & Ai for all x- In particular, there is no potential energy. In addition, 
the winding number is not quantized, and hence the defects are unstable. Finally, 
they exist only in theories with a global internal symmetry. 
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6.2. Formation 

The Kibble mechanism^^-' ensures that in theories which admit topological or 
semitopological defects, such defects will be produced during a phase transition in 
the very early Universe. 

Consider a mechanical toy model, first introduced by Mazenko, Unruh and 
Wald^^-* in the context of inflationary Universe models, which is useful in under- 
standing the scalar field evolution. Consider (see Fig. 37) a lattice of points on a 
flat table. At each point, a pencil is pivoted. It is free to rotate and oscillate. The 
tips of nearest neighbor pencils are connected with springs (to mimic the spatial 
gradient terms in the scalar field Lagrangean). Newtonian gravity creates a poten- 
tial energy V{ip) for each pencil ((/? is the angle relative to the vertical direction). 
V{ip) is minimized for \ip\ = rj (in our toy model rj = 7r/2). Hence, the Lagrangean 
of this pencil model is analogous to that of a scalar field with symmetry breaking 
potential (6.1). 



Figure 37: The pencil model: the potential energy of a simple pencil has the 
same form as that of scalar fields used for spontaneous symmetry breaking. The 
springs connecting nearest neighbor pencils give rise to contributions to the energy 
which mimic spatial gradient terms in field theory. 

At high temperatures T ^ Tc, all pencils undergo large amplitude high fre- 
quency oscillations. However, by causality, the phases of oscillation of pencils with 
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large separation s are uncorrelated. For a system in thermal equilibrium, the length 
s beyond which phases are random is the correlation length C,{t). However, since 
the system is quenched rapidly, there is a causality bound on ^: 

m<t, (6.28) 

where t is the causal horizon. 

The critical temperature Tc is the temperature at which the thermal energy 
is equal to the energy a pencil needs to jump from horizontal to vertical position. 
For T < Tc, all pencils want to lie flat on the table. However, their orientations are 
random beyond a distance of ^{t) determined by equating the free energy gained 
by symmetry breaking (a volume effect) with the gradient energy lost (a surface 
effect). As expected, C,(T) diverges at T^. Very close to Tc, the thermal energy T 
is larger than the volume energy gain Ecorr in a correlation volume. Hence, these 
domains are unstable to thermal fluctuations. As T decreases, the thermal energy 
decreases more rapidly than Ecorr- Below the Ginsburg temperature Tq, there is 
insufficient thermal energy to excite a correlation volume into the state ip = 0. 
Domains of size 

ate) ~ A-^-^ (6.29) 

freeze out ' •*. The boundaries between these domains become topological de- 
fects. 

We conclude that in a theory in which a symmetry breaking phase transitions 
satisfies the topological criteria for the existence of a fixed type of defect, a network 
of such defects will form during the phase transition and will freeze out at the 
Ginsburg temperature. The correlation length is initially given by (6.29), if the 
field (f is in thermal equilibrium before the transition. Independent of this last 
assumption, the causality bound implies that ^(tg) < to- 

For times t > Iq the evolution of the network of defects may be complicated 
(as for cosmic strings) or trivial (as for textures). In any case (see the caveats 
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of Refs. 143 and 144), the causality bound persists at late times and states that 
even at late times, the mean separation and length scale of defects is bounded by 

m < t. 

Applied to cosmic strings, the Kibble mechanism implies that at the time of 
the phase transition, a network of cosmic strings with typical step length ^(tc) will 
form. According to numerical simulations -* , about 80% of the initial energy is 
in infinite strings and 20% in closed loops. 

Note that the Kibble mechanism was discussed above in the context of a global 
symmetry breaking scenario. As pointed out in Ref. 146, the situation is more 
complicated in local theories in which gauge field can cancel spatial gradients in ip 
in the energy functional, and in which spatial gradients in ip can be gauged away. 
Nevertheless, as demonstrated numerically (in 2 + 1 dimensions) in Refs. 42 and 
147 and shown analytically in Ref. 148, the Kibble mechanism also applies to local 
symmetries. 



Figure 38: Formation of a loop by a self intersection of an infinite string. Accord- 
ing to the original cosmic string scenario, loops form with a radius R determined 
by the instantaneous coherence length of the infinite string network. 

The evolution of the cosmic string network for t > Iq is complicated (see 
Section 6.4). The key processes are loop production by intersections of infinite 
strings (see Fig. 38) and loop shrinking by gravitational radiation. These two 
processes combine to create a mechanism by which the infinite string network loses 
energy (and length as measured in comoving coordinates). It will be shown that as 
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a consequence, the correlation length of the string network is always proportional 
to its causality limit 

m-t. (6.30) 

Hence, the energy density Poo(i) in long strings is a fixed fraction of the background 
energy density pdt) 

Poo(t) ~ //{(t)"^ ~ fJ^t-'^ (6.31) 



or 

Poo(^) 
Pcit) 



Gfi. (6.32) 



We conclude that the cosmic string network approaches a "scaling solution" in 
which the statistical properties of the network are time independent if all distances 
are scaled to the horizon distance. 

Applied to textures, the Kibble mechanism implies that on all scales r > tQ, 
field configurations with winding number nw > Ucr are frozen in with a probability 
p{ncr) per volume r^. The critical winding number ricr is defined as the winding 
number above which field configurations collapse and below which they expand. 
Only collapsing configurations form clumps of energy which can accrete matter. 

The critical winding n^r was determined numerically in Refs. 149 & 150 and 
analytically in Ref. 151 (see also Ref. 152). It is slightly larger than 0.5. The 
probability p{ncr) can be determined using combinatorial arguments^^^-* . 

For t > tQ, any configuration on scale ~ t with winding number riw > ricr 
begins to collapse (before t, the Hubble damping term dominates over the spatial 
gradient forces, and the field configuration is frozen in comoving coordinates). 
After unwinding, <^(x) is homogeneous inside the horizon. 

The texture model thus also leads to a scaling solution: at all times t > tc 
there is the same probability that a texture configuration of scale t will enter the 
horizon, become dynamical and collapse with a typical time scale t. 
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6.3. Topological Defects and Cosmology 

Topological defects are regions in space with trapped energy density. By New- 
tonian gravity, these defects can act as seeds about which the matter in the Universe 
clusters, and hence they play a very important role in cosmology. 

As indicated in Table 1, theories with domain walls or with local monopoles 
are ruled out, and those with only local textures do not give rise to a structure 
formation model. As mentioned earlier, theories with domain walls are ruled out 
since a single wall stretching across the present Universe would overdose it. Local 
monopoles are alos problematic since they do not interact and come to dominate the 
energy density of the Universe. Local textures do not exist as coherent structures 
with nonvanishing gradient energy since the gauge fields can always compensate 
scalar field gradients. 

Let us demonstrate explicitly why stable domain walls are a cosmological 
disaster^^^-* . If domain walls form during a phase transition in the early Uni- 
verse, it follows by causality (see however the caveats of Refs. 143 and 144) that 
even today there will be at least one wall per Hubble volume. Assuming one wall 
per Hubble volume, the energy density pow of matter in domain walls is 

PDwit) - V^t-^ (6.33) 

whereas the critical density pc is 

Pc = H' ^ - ml,t-\ (6.34) 

Hence, for r] ~ 10^^ GeV the ratio of (6.33) and (6.34) is 

P^{t)r. ( ^]\r^t) r. 10^\ (6.35) 

Pc \mpij 

The above argument depends in an essential way on the dimension of the defect. 
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One cosmic string per Hubble volume leads to an energy density pcs in string 

Pcsr-rfr'^. (6.36) 

Later in this section we shall see that the scaling (6.36) holds in the cosmic string 
model. Hence, cosmic strings do not lead to cosmological problems. On the con- 
trary, since for GUT models with r] ~ 10^^ GeV 

P£l^ (JL\ ^10-6^ (6.37) 

Pc \mpij 

cosmic strings in such theories could provide the seed perturbations responsible for 
structure formation. 

Theories with local monopoles are ruled out on cosmological grounds^^^-* (see 
again the caveats of Refs. 143 and 144) for rather different reasons. Since there are 
no long range forces between local monopoles, their number density in comoving 
coordinates does not decrease. Since their contribution to the energy density scales 
as a~^(t), they will come to dominate the mass of the Universe, provided r] is 
sufficiently large. 

Theories with global monopoles^^"*'^^^-* are not ruled out, since there are long 
range forces between monopoles which lead to a "scaling solution" with a fixed 
number of monopoles per Hubble volume. 

In the following we will describe aspects of two of the promising topological 
defect models of structure formation, those based on cosmic strings and on global 
textures. The global monopole scenario is in many aspects similar to the texture 
theory. 
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6.4. Cosmic String Evolution and Scaling 

If the evolution of the cosmic string network were trivial in the sense that all 
strings would only stretch as the universe expands, there would be an immediate 
cosmological disaster. Consider a fixed comoving volume V with a string passing 
through. The energy in radiation decreases as a~^(t) while the energy in string 
increases as a{t). Hence trivial evolution would immediately lead to a string domi- 
nated universe, a cosmological disaster. In order to study the evolution of a cosmic 
string network, it is neccessary to know the effective action for a string, and to 
study what happens when two strings cross. 

The equations of motion of a string are determined by the Nambu action 

1/2 



S = —^ / dadr ( — det g^^^ ) a, 6 = 0, 1 



(6.38) 



(2) 

where g^^^ is the world sheet metric and a and r are the world sheet coordinates. In 
fiat space-time, r can be taken to be coordinate time, and a is an affine parameter 
along the string. In terms of the string coordinates X^(o", r) and the metric g^iJ 
of the background space-time, 

^1? = Xi:,XlgjS . (6.39) 

From general symmetry considerations, it is possible to argue that the Nambu 
action is the correct action. However, I shall follow Foerster^^^-* and Turok^^^-* and 
give a direct heuristic derivation. We start from a general quantum field theory 
Lagrangean Cqft- The action is 

s = jdScQFTiviy)) (6.40) 

We assume the existence of a linear topological defect at X^{(t,t). The idea now 
is to change variables such that a and r are two of the new coordinates, and to 
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expand S to lowest order in w/R, where w is the width of the string and R its 
curvature radius. As the other new coordinates we take the coordinates p^ and p^ 
in the normal plane to X^{cr, r). Thus the transformation takes the old coordinates 
|/^(/i = 0, 1, 2, 3) to new ones cr" = (r, a, p^, p"^): 

l/^K)=X>,r)+p^nf(a,r) (6.41) 

where i = 2,3 and nf are the basis vectors in the normal plane to the string world 
sheet. The measure transforms as 



(fy= I dadTdp^dp^{detM^) (6.42) 



with 

The determinant can easily be evaluated using the following trick 

det M^ = (- det r]^^M^M^f^ = y^-det Dab (6.44) 



Hence 



/I /9 r 
dadr (- det g^]^) J dpV^vi^, r, p\ P^)) + O (^ 

= -p / dadr (- det g^^A +o( 



(W 

R 



(6.46) 



— p is the integral of C in the normal plane of X. To first order in w/R, it equals 
the integral of — 7i; hence it is the mass per unit length. 
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This derivation of the Nambu action is instructive as it indicates a method for 
calculating corrections to the equations of motion of the string when extra fields are 
present, e.g. for superconducting cosmic strings. It also gives a way of calculating 
the finite thickness corrections to the equations of motion which will be important 
at cusps (see below). 

In flat space-time we can consistently choose t = t,x ■ x' = and x'^ + x' = 0. 
The equations of motion derived from the Nambu action then become 

X - x" = . (6.47) 

where / indicates the derivative with respect to a. The general solution can be 
decomposed into a left moving and a right moving mode -* 

x(t, a) = - [a{a -t)+ b(a + t)] (6.48) 



The gauge conditions imply 

a^ = b" = 1 (6.49) 



a2 ^-2 



For a loop, x{a, t) is periodic and hence the time average of a and b vanish, a and 
b are hence closed curves on the unit sphere with vanishing average. Two such 
curves generically intersect if they are continuous. An intersection corresponds to 
a point with x' = and x = 1. Such a point moving at the speed of light is called a 
cusp. x{a, t) need not be continuous. Points of discontinuity are called kinks. Note 
that both cusps and kinks will be smoothed out by finite thickness effects^^^-*. 

The Nambu action does not describe what happens when two strings hit. This 
process has been studied numerically for both global^^^-* and local^^^-* strings. The 
authors of these papers set up scalar field configurations corresponding to two 
strings approaching one another and evolve the complete classical scalar field equa- 
tions. The result of the analysis is that strings do not cross but exchange ends, 
provided the relative velocity is smaller than 0.9. Thus, by self intersecting, an 
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infinite string will split off a loop (Fig. 38). An important open problem is to 
understand this process analytically. For a special value of the coupling constant 
Ruback^^^-* has given a mathematical explanation (see also Shellard and Ruback 
in Ref. 161). 

There are two parts to the nontrivial evolution of the cosmic string network. 
Firstly, loops are produced by self intersections of infinite strings. Loops oscillate 
due to the tension and slowly decay by emitting gravitational radiation. Combining 
the two steps we have a process by which energy is transferred from the cosmic 
string network to radiation. 

There are analytical indications that a stable "scaling solution" (already de- 
scribed in Section ) for the cosmic string network exists. In the scaling solution, on 
the order of 1 infinite string segment crosses every Hubble volume. The correlation 
length ^(t) of an infinite string is thus of the order t. A heuristic argument for the 
scaling solution is due to Vilenkin^-*. Take z/(t) to be the mean number of infinite 
string segments per Hubble volume. Then the energy density in infinite strings is 

Poo(t)=/i/>(t)t-2 (6.50) 

The number of loops n(t) produced per unit volume is proportional to the square 
of u, since it takes two string segments to generate a string intersection. Hence, 

^ = c9H~' (6.51) 

dt ^ ' 

where c is a constant of the order 1. Conservation of energy in strings gives 

^^ + 1 p^it) = -c'fit ^ = -c'fiuH-' (6.52) 

dt 2t^^ ' ^ dt ^ ^ ' 

or, written as an equation for z/(t) 

z> - — = -cc'u^t^^ (6.53) 

Thus a u ^ 1 then z> < while if i/ <^ 1 then z> > 0. Hence there will be a stable 
solution with z/ ~ 1. 
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The precise value of v must be determined in numerical simulations. These 
simulations are rather difficult because of the large dynamic range required and due 
to singularities which arise in the evolution equations near cusps. In the radiation 
dominated epoch, v is still uncertain by a factor of about 10. The ffist results 
were reported in Ref. 163. More recent results are due three groups. Bennett and 
Bouchet^^^-* and Allen and Shellard^^^-* are converging on a value 10 < i/ < 20, 
whereas Albrecht and Turok^^^-* obtain a value which is about 100. 

The scaling solution for the infinite strings implies that the network of strings 
looks the same at all times when scaled to the Hubble radius. This should also im- 
ply that the distribution of cosmic string loops is scale invariant in the same sense. 
At present, however, there is no convincing evidence from numerical simulations 
that this is really the case. 

A scaling solution for loops implies that the distribution of Ri{t), the radius of 
loops at the time of formation, is time independent after dividing by t. To simplify 
the discussion, I shall assume that the distribution is monochromatic, i.e. 

Ri{t)/t = a. (6.54) 

Based on Fig. 38, we expect a ~ 1. The numerical simulations^^^"^^^-*, however, 
now give a < 10^^. 

From the scaling solution (6.50) for the infinite strings we can derive the scaling 
solution for loops. We assume that the energy density in long strings - inasmuch 
as it is not redshifted - must go into loops. [3 shall be a measure for the mean 
length £ in a loop of "radius" R 

e = (3R. (6.55) 

If per expansion time and Hubble volume about 1 loop of radius Ri{t) is produced, 
then we know that the number density in physical coordinates of loops of radius 
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Ri{t) is 

n{Ri{t),t) = ct-^ (6.56) 

with a constant c which can be calculated from (6.50), (6.54) and (6.55). Neglecting 
gravitational radiation, this number density simply redshifts 

n{R,t)= (j^^^) n(i?,t;(i?)), (6.57) 

where tf{R) is the time when loops of radius R are formed. Isolating the R 
dependence, we obtain 

n{R,t) ^ R'^z{R)-^ (6.58) 

where z{R) is the redshift at time t = R. We have the following special cases: 

n(i?, t) ~ i?"^/ V3/2 t<teq 

niR, t) ~ R-^'hl'q^t-'^ t>teq, tf{R) < Uq (6.59) 

n{R, t) ~ i?"2^-2 t > teg , tf{R) > teq ■ 

The proportionality constant c is 

c=-p-^a-^i) (6.60) 

(see e.g. Ref. 167). In deriving (6.60) it is important to note that n{Ri{t),t)dRi 
is the number density of loops in the radius interval [Ri, Ri + dRi]. Hence, in the 
radiation dominated epoch 

with 

u = ^(3~^a^/'^9. (6.62) 

From (6.62) we can read off the uncertainties in z/ based on the uncertainties 
in the numerical results. Both a ' and i> are determined only up to one order of 
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magnitude. Hence, any quantitative results which depend on the exact value of v 
must be treated with a grain of salt. 

Gravitational radiation leads to a lower cutoff in n{R,t). Loops with radius 
smaller than this cutoff were all formed at essentially the same time and hence have 
the same number density. Thus, n{R) becomes flat. The power in gravitational 
radiation Pq can be estimated using the quadrupole formula^^®-' . For a loop of 
radius R and mass M 

Pg = ^G<QQ>, (6.63) 

where Q is the quadrupole moment, Q ~ MR^^ and since the frequency of oscilla- 
tion is cj = R~^ 

Pg ~ G{MR^fu^ ~ (G/i)/i . (6.64) 

Even though the quadrupole approximation breaks down since the loops move 
relativistically, (6.64) gives a good order of magnitude of the power of gravitational 
radiation. Improved calculations give^^^-* 

Pg = liGfi)fi (6.65) 

with 7 ~ 50. (6.55) and (6.65) imply that 

R = ^Gfj. (6.66) 

with 7 = 7//3 ~ 5 (using P ~ 10). Note that the rate of decrease is constant. 
Hence, 

R(t) = Ri-{t- ti)^Gfx (6.67) 

and the cutoff loop radius is 

Re ~ ^Gfiti . (6.68) 

Let us briefly summarize the scaling solution 
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1) At all times the network of infinite strings looks the same when scaled by the 
Hubble radius. A small number of infinite string segments cross each Hubble 
volume and pooit) is given by (6.50). 

2) There is a distribution of loops of all sizes < R < t. Assuming scaling for 
loops, then 



-4 / ^(t) 



n{R, t) = z/i?"* ( ^^ 1 , Re [^Gfit, at] (6.69) 

where a~^R is the time of formation of a loop of radius R. Also 

n{R, t) = n{^Gfit, t) , R < ^Gfit . (6.70) 

Although the qualitative characteristics of the cosmic string scaling solution 
are well established, the quantitative details are not. The main reason for this 
is the fact that the Nambu action breaks down at kinks and cusps. However, 
kinks and cusps inevitably form and are responsible for the small scale structure 
on strings. In fact, coarse graining by integrating out the small scale structure 
may give an equation of state for strings which deviates from that of a Nambu 
string^™-*. Attempts at understanding the small scale structure on strings are at 
present under way^"^^-*. 

6.5. Cosmic Strings and Structure Formation 

The starting point of the structure formation scenario in the cosmic string 
theory is the scaling solution for the cosmic string network, according to which at 
all times t (in particular at tf,q, the time when perturbations can start to grow) 
there will be a few long strings crossing each Hubble volume, plus a distribution 
of loops of radius R <^t (see Fig. 39). 

The cosmic string model admits three mechanisms for structure formation: 
loops, filaments, and wakes. Cosmic string loops have the same time averaged field 
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Figure 39. Sketch of the scahng solution for the cosmic string network. The box 
corresponds to one Hubble volume at arbitrary time t. 

as a point source with mass"*^ ' 

M{R)=(3Rid, (6.71) 

R being the loop radius and j3 ~ 2n. Hence, loops will be seeds for spherical 
accretion of dust and radiation. 

For loops with R <tf,q, growth of perturbations in a model dominated by cold 
dark matter starts at tgg. Hence, the mass at the present time will be 

M{R,to) = z{teq)(3Ri2. (6.72) 

In the original cosmic string model^^'^^'^^-* it was assumed that loops dominate 
over wakes. In this case, the theory could be normalized {i.e., /i could be deter- 
mined) by demanding that loops with the mean separation of clusters dd (from 
the discussion in Section 6.4 it follows that the loop radius R{dci) is determined 
by the mean separation) accrete the correct mass, i.e., that 

M{R{dd)M) = l^^^MQ. (6.73) 

109 



This condition yields^^) 

/i ^ lO^^GeV^ (6.74) 

Thus, if cosmic strings are to be relevant for structure formation, they must arise 
due to a symmetry breaking at energy scale r] ~ 10 GeV. This scale happens to 
be the scale of unification (GUT) of weak, strong and electromagnetic interactions. 
It is tantalizing to speculate that cosmology is telling us that there indeed was new 
physics at the GUT scale. 



Figure 40. Sketch of the mechanism by which a long straight cosmic string S 
moving with velocity v in transverse direction through a plasma induces a velocity 
perturbation Av towards the wake. Shown on the left is the deficit angle, in the 
center is a sketch of the string moving in the plasma, and on the right is the sketch 
of how the plasma moves in the frame in which the string is at rest. 

The second mechanism involves long strings moving with relativistic speed in 
their normal plane which give rise to velocity perturbations in their wake^^"^-*. The 
mechanism is illustrated in Fig. 40: space normal to the string is a cone with 
deficit angle ^^^-^ 

a = 8iTGfi. (6.75) 

If the string is moving with normal velocity v through a bath of dark matter, a 
velocity perturbation 

6v = AnGfj.v)' (6.76) 

[with 7 = (1 — f^)~^'^] towards the plane behind the string results. At times af- 
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ter teq, this induces planar overdensities, the most prominent {i.e., thickest at the 
present time) and numerous of which were created at teq, the time of equal mat- 
ter and radiation^^'^^'^^-* . The corresponding planar dimensions are (in comoving 
coordinates) 

teqZ{teq) X teqZ{teq)v ~ (40 X AOv) MpC^ . (6.77) 

The thickness d of these wakes can be calculated using the Zel'dovich approximation -*.| 
The result is 

d ~ Gl2V-f{v)z{teqf teq ~ Av MpC . (6.78) 

Wakes arise if there is little small scale structure on the string. In this case, 
the string tension equals the mass density, the string moves at relativistic speeds, 
and there is no local gravitational attraction towards the string. 

In contrast, if there is small scale structure on strings, then the string tension 
T is smaller^™-* than the mass per unit length /i and the metric of a string in z 
direction becomes^^^-* 

ds"^ = (1 + hQo){dt'^ - dz^ - dr^ - (1 - 8Gfi)r^dy^) (6.79) 

with 

hoo = AGifi - T) In — , (6.80) 

tq being the string width. Since /iqo does not vanish, there is a gravitational 
force towards the string which gives rise to cylindrical accretion, thus producing 
filaments. 

As is evident from the last term in the metric (6.79), space perpendicular to 
the string remains conical, with deficit angle given by (6.75). However, since the 
string is no longer relativistic, the transverse velocities v of the string network are 
expected to be smaller, and hence the induced wakes will be shorter and thinner. 
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Which of the mechanisms - filaments or wakes - dominates is determined by 
the competition between the velocity induced by /iqo and the velocity perturbation 
of the wake. The total velocity is^^^-* 



u = _ 2^^(/^ ^) _ AnGiivMv) , (6.81) 

f7(f) 



the first term giving filaments, the second producing wakes. Hence, for small v the 
former will dominate, for large v the latter. 

By the same argument as for wakes, the most numerous and prominent fila- 
ments will have the distinguished scale 

teqZ{teq) X df X df (6.82) 

where df can be calculated using the Zel'dovich approximation. 

The cosmic string model predicts a scale-invariant spectrum of density pertur- 
bations, exactly like inflationary Universe models but for a rather different reason. 
Consider the r.vn.s. mass fluctuations on a length scale 2'Kk^^ at the time tH{k) 
when this scale enters the Hubble radius. From the cosmic string scaling solution 
it follows that a fixed (z.e., tuik) independent) number v of strings of length of 
the order t}j{k) contribute to the mass excess 6M{k, t}j{k)). Thus 

^ (A:, tnm ^ ^^l^^^l ^vG^. (6.83) 

Note that the above argument predicting a scale invariant spectrum will hold for 
all topological defect models which have a scaling solution, in particular also for 
global monopoles and textures. 

The amplitude of the r.m.s. mass fluctuations (equivalently: of the power 
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spectrum) can be used to normalize Gjj,. Since today on galaxy cluster scales 

{k, to) ~ 1 , (6.84) 



M 

the growth rate of fluctuations linear in a{t) yields 

SM 



M 
and therefore, using v ~ 10, 



{k,te,) ^lO-\ (6.85) 



-5 



G/i~10"\ (6.86) 



A big advantage of the cosmic string model over inflationary Universe models 
is that HDM is a viable dark matter candidate. Cosmic string loops survive free 
streaming, as discussed in Section 3.4, and can generate nonlinear structures on 
galactic scales, as discussed in detail in Refs. 61 and 62. Accretion of hot dark 
matter by a string wake was studied in Ref. 63. In this case, nonlinear perturba- 
tions develop only late. At some time t^i, all scales up to a distance qmax from the 
wake center go nonlinear. Here 

gmax ~ G^lV'-f{v)z{teq)heq ~ 4w MpC , (6.87) 

and it is the comoving thickness of the wake at tni- Demanding that tni corresponds 
to a redshift greater than 1 leads to the constraint 

G'/i>5■10"^ (6.88) 

The power spectra in the cosmic string models with CDM and HDM are obvi- 
ously different on scales smaller than the maximal neutrino free streaming length 
(3.27). Recent calculations^^^'^^^-* of the power spectra are shown in Fig. 41. 
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Figure 41: Power spectra for cosmic string HDM and CDM theories (dashed 
curves), compared to those for inflationary HDM and CDM models (solid curves). 
In each case, the top curve is for CDM, the bottom one for HDM. Note that there 
is substantial power on small scales in the cosmic string HDM theory. 

6.6. Global Textures and Structure Formation 

The starting point of the texture scenario of structure formation^^-* is the scal- 
ing solution for textures: at any time t, there is a fixed probability p{nyj) driw 
that the scalar field configuration over a Hubble volume covers between n^ and 
flu, + dn^i, of the vacuum manifold, i.e., we have a texture with winding number in 
the interval [n^o, n^ + dn^] entering the Hubble radius. 

The dynamics of a texture is easy to understand. Consider the spherically 
symmetric texture configuration of (6.24) with x(r) increasing from to Xmax over 
a distance d. If d is larger than the Hubble radius, then the Hubble damping 
term dominates the equation of motion for ip and the field configuration is frozen 
in. Once the Hubble radius t catches up with d, the microphysical forces become 
dominant and the texture field begins to evolve. 
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Figure 42: A sketch of the forces acting on a spherically symmetric texture con- 
figuration and which cause unwinding in case (a) in which the winding number is 
larger than the critical winding, and dissipation if the winding is smaller than its 
critical value (case (b)). r is the distance from the center of the texture, and the 
vertical axis shows the value of the x field. 

The evolution of (/? tends to minimize the field energy. Consider first large 
distances from the texture center. The spatial gradient energy can be decreased 
by having Xmax increase (if Xmax > Xc) or decrease (if Xmax < Xc) (see Figure 42). 
The winding associated with Xc is called the critical winding ric (see (6.26)). For 
a single texture in an infinite volume we would expect 

Xc = |(^.e.,nc = 0.5). (6.89) 

For realistic textures there will be a "finite volume cutoff' determined by the sepa- 
ration of textures. A semi-analytical analysis and numerical simulations give^^^"^^^-* 

nc~ 0.6. (6.90) 

If nw < ric, then the field configuration will relax to a trivial one. No local- 
ized energy concentrations will be generated, and we cannot speak of a "texture." 
However, if nw > ric the field evolution will be more interesting. At large r, x(r) 
will increase. In addition, the radius r{x) where x takes on a fixed value x tends 
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to decrease, since this leads to a concentration of gradient energies over a smaller 
region. Hence, the field configuration will contract (see Fig. 42), with increasing 
total winding number. Eventually, close to r = there is sufficient tension energy 
for ip to be able to leave the vacuum manifold and jump from x = to x = ^• 
This is the texture unwinding event. After unwinding, energy is radiated radially 
in the form of Goldstone bosons. 

In the texture model it is the contraction of the field configuration which leads 
to density perturbations^'^®-'. At the time when the texture enters the horizon, an 
isocurvature perturbation is established: the energy density in the scalar field is 
compensated by a deficit in radiation. However, the contraction of the scalar field 
configuration leads to a clumping of gradient and kinetic energy at the center of 
the texture (Fig. 43). This, in turn, provides the seed perturbations which cause 
dark matter and radiation to collapse in a spherical manner^^^'^®^-*. 



Figure 43: A sketch of the density perturbation produced by a collapsing texture. 
The left graph shows the time evolution of the field x(r) as a function of radius r 
and time (see (5.18)). The contraction of x(^) leads to a spatial gradient energy 
perturbation at the center of the texture, as illustrated on the right. The energy is 
denoted by p. Solid lines denote the initial time, dashed lines are at time t + At, 
and dotted lines correspond to time t + 2At, where At is a fraction of the Hubble 
expansion time (the typical time scale for the dynamics). 
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As in the cosmic string model, also in the global texture scenario the length 
scale of the dominant structures is the comoving Hubble radius at tgg- Textures 
generated at tgq are the most numerous, and the perturbations induced by them 
have the most time to grow. 

As mentioned in the previous subsection, the texture model predicts a scale- 
invariant spectrum of density perturbations. Hence, in order to differentiate topo- 
logical defect models from inflationary scenarios, and to distinguish between differ- 
ent topological defect theories, we need statistics which are not determined by the 
power spectrum alone. We need statistics which are sensitive to the non-random 
phases of topological defect models. One such statistic is the genus curve^^^-* . For 
a surface S embedded in R^, the genus g is 

g{S) = # of holes oi S — # of disconnected components of 5 -|- 1 . (6.91) 

The genus g can now be evaluated for the isodensity surface S{p), the surface of 
points in space with density equal to p. The curve 

g{p) = g{S{p)) (6.92) 

is the genus curve. To reduce numerical errors, g can also be evaluated based on a 
cell decomposition of the volume. Now, g{n) is the genus of the boundary of the 
cell complex in which each cell contains more than n galaxies. In this case, the 
genus is simply 

g = l-^(V-E-F) (6.93) 

where V, E, F are the number of vertices, edges, and faces of the polygonal surface, 
respectively. 

As shown in Fig. 44, the genus statistic is able to distinguish between models 
with the same power spectrum but different phase correlations and topology^^^-*. 
For a texture toy model, the genus curve is mostly negative, for a cosmic string 
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Figure 44: Comparison of tlie genus curve (genus as a function of galaxy density) 
of different toy models of structure formation. Except for the Gaussian model, all 
theories have the same linear power spectrum. The 'filament', 'wake' and 'texture' 
toy models are based on laying down at random linear, planar and spherical over- 
dense regions of galaxies. Thus, the figure demonstrates that the genus statistic is 
able to distinguish between theories with different topologies but identical power 
spectra. The 'CDM' model predictions are computed from linear theory, and the 
'Poisson' model is obtained by randomly distributing galaxies. See the senior thesis 
by Aguirre^^^-* for further details. 

wake model, it is predominately positive. The differences compared to a random 
phase inflationary model are statistically significant. 

The differences shown in Fig. 44 will only be apparent in large-scale samples of 
galaxies, i.e., on scales exceeding the comoving radius at teq. Such samples should, 
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however, become available in the near future, and at that point genus curve and 
other statistics sensitive to non-random phases should become a powerful tool for 
distinguishing the predictions of the different models of structure formation. 

A final word concerning textures: since they are short-lived, only CDM is a 
viable dark matter candidate in the context of this structure formation scenario. 

7. Cosmic Microwave Background Anisotropics 

As mentioned in Section 3, the near-isotropy of the CMB is the strongest evi- 
dence in support of the cosmological principle. By the same reasoning, any density 
inhomogeneities in the early Universe will give rise to CMB anisotropies. Since 
our present theories of galaxy formation are based on the gravitational instability 
scenario, they predict such inhomogeneities. The CMB temperature fiuctuations 
probe the structure of space at tree, the time of last scattering, a time when the 
density perturbations still have a small amplitude and can be analyzed in linear 
theory. Hence, a study of CMB anisotropies will yield a lot of constraints for 
structure formation models. The information gained will be robust, i.e., indepen- 
dent of the uncertainties of nonlinear gravitational and hydrodynamical effects, 
but it will deal only with large scales (comparable or larger than the comoving 
horizon at tree)- In this section we shall give a brief overview of the theory of CMB 
anisotropies and summarize some recent observational results. 

7.1. Basics 

As illustrated in Fig. 45, there are three main sources of CMB anisotropies. 
The first are gravitational potential perturbations at t^ee which lead to fiuctuations 
of the surface of last scattering. This produces deviations in the light travel time 
between last scattering and detection, and - given that the photons have the same 
temperature on the surface of last scattering - to temperature fiuctuations for the 
observer. 
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Figure 45: Space-time plot sketcliing tlie origin of CMB temperature anisotropies. 
The surface labelled Tree is the last scattering surface. O is the observer at the 
present time measuring photons 7 impinging from directions in the sky separated 
by an angle 6. The shaded area labelled C is the world volume of a local overdensity, 
leading to distortions of geodesies. Possible velocities of observer and emitter are 
indicated by Vq and Vg, respectively. 

The second source is due to gravitational perturbations along the line of sight 
which lead to deviations of the geodesies and hence to temperature differences. 
A Newtonian way of understanding this effect is to consider a photon passing 
through a large mass concentration. On the way towards the center, the photon is 
falling into a potential well and acquires a blueshift, whereas on its way out it is 
redshifted. In an expanding background, this redshift does not exactly cancel the 
initial blueshift, and a temperature fluctuation results. 

The third source contributing to CMB anisotropies are peculiar velocities on 
the surface of last scattering and of the observer. The peculiar motion of the earth 
gives rise to a dipole anisotropy^^^-* 

rry-i 

T^ I dipole V • / 

Peculiar velocities induce temperature fluctuations by means of the Doppler effect. 
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For linear adiabatic density perturbations in a matter dominated Universe, the 
line of sight contributions to 5T /T can be written as total time derivative and thus 
reduces to a contribution from the surface of last scattering and can be simply 
combined with the potential fluctuations at tree- This is the case first studied by 
Sachs and Wolfe, and the combined effect is now called the Sachs- Wolfe effect ^^'^•'. 

An analysis of the Sachs- Wolfe effect reveals a very simple relationship be- 
tween temperature fluctuations 5T /T {d) on an angular scale d and the magnitude 
of density perturbations on the corresponding lengths scale A('i9), where at last 
scattering A(i9) equals the distance subtended by two light rays with angular sepa- 
ration tD (see Fig. 45). A simple derivation^^'^^^-* of this relationship makes use of 
the gauge invariant theory of cosmological perturbations described in Section 4.4. 

The starting point is the phase space distribution function /(a;", pi) which 
would be a function of p/T exclusively in the absence of inhomogeneities. In the 
presence of inhomogeneities, the deviation of / from homogeneity is associated 
with temperature fluctuations: 

/(a:", Pi) = Kp/T + 6T) , (7.2) 

where T is the average temperature and f{p/T) is the background phase space 
density. 

The phase space distribution function satisfies the coUisionless Boltzmann 
equation 

^1^ + ?^I^ = (7-3) 

dr] ox" dr] opi 

where, as in Section 4.4, the variable r] denotes conformal time. This equation can 
be integrated along the perturbed geodesies which are given by 

— = 2 ^ (7 4) 
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and 



with 



^ = r(l + 2$), (7.5) 



r = --Pi (7.6) 

p 

and 

P^ = PiPi ■ (7.7) 

Inserting these relations into the Boltzman equation gives 

l^ + rs.) ^ = -2ra.$. (7.8) 

Since in the matter dominated period (9,^$ = we can rewrite (7.8) as 

| + /A(f + 2*)^0. (7^9) 



which imphes that 



along the perturbed geodesies. 



rry-i 

— + 2$ = const (7.10) 



For isothermal primordial perturbations (^ (tree) = 0), the result (7.10) im- 
plies that 

— (r^o) = 2^{r]rec) + l"vi{r]rec) (7.11) 

whereas for primordial adiabatic perturbations (vanishing initial entropy pertur- 
bations) 

^(^O) = -HVrec) + l'viir]rec) , (7.12) 

i.e., the combination of initial curvature fluctuations and line of sight effects leads 
to a partial cancellation of the anisotropy. The second term on the r.h.s. of (7.11) 
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and (7.12) is the Doppler term, and it arises from a determination of the constant 
in (7.10) based on considering the initial conditions at tree- 
Since $ is constant both between teq and tree and while outside the Hubble 
radius, and since 

HtH) r^ ^-^ (tn) (7.13) 

P 

at Hubble radius crossing tn, our results imply that (modulo Doppler terms) for 

adiabatic perturbations 



^(^,to) = ^$(A(^),W~^^ 



(^, to) = ^HH^), teg) ~ 77 — (A(^), tn)) . (7.14) 



We conclude that the spectrum of primordial mass perturbations can be normalized 
by CMB anisotropy detections. For a scale invariant spectrum of density pertur- 
bations, the r.m.s. temperature fluctuations are predicted to be independent oi i) 
on angular scales larger than the Hubble radius at tree (between 1 and 2 degrees). 

7.2. Specific Signatures 

All theories of structure formation give rise to Sachs- Wolfe type temperature 
fluctuations given by (7.13) and (7.14). In topological defect models there are, 
in addition, speciflc signatures which cannot be described in a linear perturbative 
analysis. 

As described in Section 6.5, space perpendicular to a long straight cosmic 
string is conical with deflcit angle given by (6.75). Consider now CMB radiation 
approaching an observer in a direction normal to the plane spanned by the string 
and its velocity vector (see Fig. 46). Photons arriving at the observer having 
passed on different sides of the string will obtain a relative Doppler shift which 
translates into a temperature discontinuity of amplitude^^^-* 

— = A7rGi^v^{v), (7.15) 

where v is the velocity of the string. Thus, the distinctive signature for cosmic 
strings in the microwave sky are line discontinuities in T of the above magnitude. 
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Figure 46: Sketch of the Kaiser-Stebbins effect by which cosmic strings produce 
hnear discontinuities in the CMB. Photons 7 passing on different sides of a moving 
string S (velocity v) towards the observer O receive a relative Doppler shift due 
to the conical nature of space perpendicular to the string (deficit angle a). 

Given ideal maps of the CMB sky it would be easy to detect strings. However, 
real experiments have finite beam width. Taking into account averaging over a 
scale corresponding to the beam width will smear out the discontinuity, and it 
turns out to be surprisingly hard to distinguish the predictions of the cosmic string 
model from that of inflation-based theories using quantitative statistics which are 
easy to evaluate analytically, such as the kurtosis of the spatial gradient map of 
the CMRis^). 

Textures produce a distribution of hot and cold spots on the CMB sky with 
typical size of several degrees •*. This signature is much easier to see in CMB 
maps. The mechanism which produces these hot and cold spots in the CMB is 
illustrated in Fig. 47. 

Photons arriving at the observer having passed through a texture as in the 
case of the ray 71 in Fig. 47 will be redshifted relative to the average photons since 
they have to climb out of a potential well, whereas those in orientation 72 will be 
blueshifted since they fall into a potential well. Taking into account reionization 
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Figure 47: Space-time diagram of a collapsing texture. The unwinding occurs at 
the point TX. The shaded areas correspond to overdense regions. Photons hke 71 
are redshifted, those hke 72 are blueshifted. 

produced by texture coUapse gives an amphtude of 5T/T of^^^'^^'^'' 

— ~ 0.06 X leTrG'?]^ . (7.16) 

A number of about ten hot and cold spots of angular scale 10° is predicted by the 
texture model. 

Theories of structure formation can now be normalized from CMB anisotropy 
data and from large-scale structure considerations. An inflationary model with 
CDM yields agreement between these two normalizations provided^^^-* 

6-1, (7.17) 

where b is the bias factor determining the ratio of fractional mass to light pertur- 
bations on a scale of 8h~^ Mpc. 



6L 

T 



, SM 
b 



(7.18) 

8h-iMpc 



8h-iMpc M 

However, agreement between galaxy and cluster correlation properties seem to 
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require^^^) 



6~2. (7.19) 



Normalizations of the texture model from large-scale structure and CMB observations^^^'^^'^'^^^-'l 
require a bias 

6-3, (7.20) 

whereas for cosmic strings the two normalizations agree well. Based both on numer- 
ical simulations and analytical calculations, a normalization of the cosmic string 
model from the COBE CMB anisotropy data gives^^^'^^^-* 

G/i= (1.3±0.5)10"^ (7.21) 



7.3. Experimental Results 

Over the past couple of years there has been a spectacular breakthrough on 
the observational front. The DMR experiment on the COBE satellite •* has pro- 
duced a temperature map of the entire sky with beam width of 7°, which shows 
a clear detection of CMB anisotropics. Independent confirmation has come from 
two 5° experiments, FIRAS^^^-* which has mapped 1/4 of the sky, and the Tenerife 
experiment ^^^-^ which surveyed a strip of 70° length in right ascension at a dec- 
lination 40°. The FIRAS data cross correlate very well with the COBE results, 
and there is even good agreement in the location of a pronounced feature in the 
Tenerife map with a that of a comparable feature in the two-year COBE maps. 

In addition, there are many small angular scale experiments which have de- 
tected anisotropics. A partial list of observational results is given in Table 2. In 
this table, "Angular Scale" denotes the beam width, the "results for bT jT'' stands 
for the variance of bT computed from the CMB maps, "cover" indicates the area 
of the sky mapped. MAX 1 and MAX 2 denote two separate MAX measurements 
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of 5T /T ^ one in a region of the sky /x Peg, the second near GUM. OVRO 1 is the 
first Owens Valley experiment, a measurement near the North Galactic Cap, the 
second is a ring survey. The large numer of anisotropy experiments which have 
announced detections of temperature fluctuations since April 1992 indicates the 
rapid progress in this field. 

To a first approximation, the present experimental results are in agreement 
with the predictions of a scale invariant spectrum of density perturbations. A 
popular way to show the results is to expand T(n) in spherical harmonics 



T(n) = Y.Y. ^imYU^). (7.22) 

I m——l 

where n is a unit vector on the sky, and to calculate the temperature correlation 
function 

< T(ni)T(n2) >= -^ J](2/ + l)QPKni " m) (7.23) 



Air 
I 



where 



< a*^ai>m' >= Ci6ii'6mm' ■ (7.24) 

For a power spectrum of density perturbations 

P{k) ~ A;" (7.25) 

the prediction for the Sachs- Wolfe contribution to ST is 

fCi ~ l"-^ (7.26) 

on scales larger than the Hubble radius at tree {i-e., for small values of /). 

A direct comparison between theory and experiment is complicated by two 
effects: the Doppler contribution to ST/T creates a peak in the I'^Ci curve at 
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values of / which correspond to wavelengths comparable to the Hubble radius at 
tree, whose amplitude depends strongly on the ionization history of the Universe. 
Reionization also leads to a decrease in C/ for large /. 

The COBE results combined with Tenerife observations favor^^®-* a value of n 
larger than what is predicted by simple inflationary models. However, the error 
bars are large and the difference is not (yet) statistically significant. At present 
there is the intriguing puzzle as to why the signal of certain small scale experiments 
is larger than the upper limit of other observations at the same angular scale 
elsewhere in the sky. A search for possible non-Gaussian features in the CMB sky 
will have high priority in the next years. 

8. Conclusions 

Modern cosmology has led to the development of several theories of structure 
formation, most prominently theories based on inflation, and topological defect 
models. These new theories are all based on the union between particle physics 
and general relativity. The models of structure formation obey the usual causality 
principle of relativistic physics. 

All of the current theories of structure formation have their problems. Most 
importantly, they do not address the cosmological constant problem but rather, 
inasmuch as they make use of scalar matter fields, make the problem worse. The 
inflationary Universe scenario is still lacking a convincing realization. Present 
versions require very special scalar field potentials. Topological defect models, 
on the other hand, do not explain why the Universe is nearly homogeneous and 
spatially fiat. In my opinion, we should regard our current theories as toy models 
with which we work and from which we learn, but which will eventually be replaced 
by improved and more convincing theories. 

Nevertheless, our present theories are predictive. To a first approximation, 
they all predict a scale invariant spectrum of density perturbations and induced 
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TABLE 2: CMB Anisotropy Results 
Experiment Angular Scale Result for -jr Cover 



Location 



COBE-DMRiy* 


3)7° 


1.1 ± 0.2 


47r 


space 


Tenerifei98) 


5.6° 


1.7 ± 0.4 


350 deg2 


ground 


FIRS197) 


4° 


1-3 


71 


balloon 


SK93199) 


1.45° 


1.4 ± 0.5 




ground 


SP91200) 


1.4° 


1.1 ± 0.5 


13.8 deg2 


ground 


ARGO201) 


1° 


2.2 ± 0.8 




balloon 


Pytlion202) 


0.75° 


3 


8deg2 


ground 


MAX 1203) 


0.5° 


<3 




balloon 


MAX 2204) 


0.5° 


4.9 ± 0.8 


6deg2 


balloon 


MSAM205) 


0.47° 


1.6 ± 0.4 


6deg2 


balloon 


White Dish206) 


0.2° 


<2.3 




ground 


OVRO 1207) 


1.8' 


< 1.9 


0.03 deg2 


ground 


OVRO 2208) 


1.8' 


3.4 ± 1.1 


0.1 deg2 


ground 



CMB anisotropies. Typically, the models contain one intrinsically free parameter 
(plus maybe a couple more parameters with which we can describe our ignorance of 
the detailed evolution of the models). The free parameter can be normalized from 
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any one of several observables. It is remarkable that the different normalizations 
of the models are consistent (to a first approximation). This lets us entertain the 
hope that we are on the right track: structure formation proceeds via gravitational 
instability (Sections 4.3 and 4.4) with the seed perturbations being provided by a 
particle physics theory of the very early Universe. 

There is already a wealth of observational data which is fit quite well by our 
present toy models. More and higher accuracy data is rapidly becoming available. 
The data concerns on one hand structure in the Universe gleamed from optical and 
infrared galaxy surveys, and on the other hand from the temperature map of the 
CMB sky. 

With the wealth of data available and steady flow of new observational results, 
and given that many important questions remain unresolved, modern cosmology 
will remain an exciting area of research for the forseeable future. 
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